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ABSTRACT 


A simplified scheme Is presented, for the heating of clear and cloudy 
air by solar and infrared radiation transfer, designed for use in tropospheric 
general circulation models with coarse vertical resolution. 

A bulk transmission function is defined for the infrared transfer. The 
interpolation factors, required for computing the bulk transmission function, 
are parameterized as functions of such physical parameters as the thickness 
of the layer, the pressure, and the mixing ratio at a reference level . 

The computation procedure for solar radiation is significantly simpli- 
fied by the introduction of two basic concepts. The first is that the solar 
radiation spectrum can be divided into a scattered part, for which Rayleigh 
scattering is significant but absorption by water vapor is negligible, and an 

absorbed part for which absorption by water vapor is significant but Rayleigh 

| 

scattering is negligible. The second concept is that of an equivalent cloud 
water vapor amount which absorbs the same amount of radiation as the cloud. 
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(fl*3) 


where c|jj q ‘ s the mixing ratio of water vapor* A similar expression for the effective 
carbon 2 dioxide amount is given by equation (IV* 17). 


Insert 2, p, 11-2 


Following Yamamoto (1952), the total transmission function of a mixture of 
wafer vapor and carbon dioxide is assumed to be represented by the product of their 
respective transmission functions, that is 


t{u / T) - t Ho q( u h o O' T) t CO„ (u CO„' T) 


* 

9 -^9 


and 


t co5' j co„' 


» 2 ~ 1 « 2 '“' ~ 2 

As discussed in section 3(ii), t cq (v^q / T) is assumed tojse constant 
between any two isobaric levels. 2 2 Thus t and t vary only 

with T„ n and t., n , respectively. Consequently, unless required for clarl by, 
the 2 U subscript FLO will henceforth be dropped from q 

- d V 


Insert 3, p, 11-9 .. 

Experiments have been conducted with the GCM to determine, once and for 
all, the functional dependence expressed by equation (11.22). In these experiments, 
the integrals appearing in (11.16) and (IIJ7) were evaluated numerically by dividing 
the GCM layer under consideration into thin sublayers of 10 mb thickness. In this way. 



Insert 4, p . 11-15 




Following Yamamoto (1952), the total transmission function of a mixture 
of water vapor and carbon dioxide is assumed to be represented by the product of 
their respective transmission functions, that is 


and 


t(u h 2 o' u CO r t) 

%V Uc *cy T) 


T H 2 0 (u H 2 0' T) t C0 2 (u C0 2 ,1) 
= t H 2 0 (u H 2 0' T) t C0 2 (u C 0 2 ' t) 


0 


Yamamoto presents values of q, 


h 2 o 


and 


T co, 


(denoted in Tables 1 1 —2 to 


II-5 by Obs.) that are based upon experimental laboratory measurements of a 
generalized absorption coefficient in certain discrete spectral ranges. Using this 
experimental data, we deduce the following empirical equations for the transmission 
functions of water vapor and carbon dioxide. 


Insert 5, p. IV— 1 

The effective water vapor amount in an air column of height ^ , u 
can be obtained from equation (II. 3) if the continuous vertical distribution of water 
vapor mixing ratio q is known. Since the UCLA 3~layer GCM predicts q from the 
moisture conservation equation only at discrete levels, some approximation technique 
must be employed for equation (11.3), Smith (1966) has shown that the seasonal and 
fdtitudinal northern hemisphere mean value of the mixing ratio in the lower troposphere 
can be expressed by 

‘ q=q o (f-)Y / (IV. 1) 

o 

where q is a known value of the mixing ratio at some pressure p , Utilizing 
different forms of this relation. In conjunction witfi the values of 0 q predicted 
by the GCM, two methods for the evaluation of u will be developed subsequently. 
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.To accurately compute the fluxes of long wave radiation and solar ' 

V*a ^ - . - 

TfftsekriHen In the atmosphere. It Is necessary to discretize the spectrum and apply 

the equation of radiative transfer to each narrow v/avelength intervale Further- 

more, the atmosphere must be divided Into a large number of thin layers for -the: 

^computation^even when the general circulation model (GCM) has only a few 

layers. The amount of computer time required to calculate the radiative fluxes 

throughout the global atmosphere by this method is enormous and would 

exceed the^tlme required for al leather computations of the GCM-eombmeet*. 

Consequently, it is desirable to have a simplified radiation computation scheme 
* ' usc.i~£\ 

for a GCM coarse vertical resolution. 


C— *i c/doidS'ez c:c ,"c:'‘ *■ T" : - 

. The of the atmosphere strongly infiuencejf the radlaMonal 

heating field. When the sky Is covered by a dense cloud layer, the net flux of 

long wave radiation in the cloud layer Is zero^-TFr^£op and bottom of the cloud 

layer are trea<t©d=as black body radiators at their respective temperatures. The 

cooling due to long wave radiation Is therefore significantly modified throughout 

1><? <2 

the atmospheric column containing the cloud layer. Furthermore, sub- 

A 

sfantiai fraction of the Incident solar radiation is reflected or absorbed by a cloud 

layer, the vertical distribution of heating solar radiation Is taraeby affected 

A A A 

by the cloud layer. 
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Ho LONG WAVE RADIATION 


1 » Basic equations 

The atmospheric long wave (thermal) radiation with wavelengths 

between 2.5 p and 40 p Is m ainly affe c ted by the absorption by water vapor, 

. . Ab ' - a /7 ey\ " . . ... ~ - 

carbon dioxide and ozone ^catfering A by air molecules is negligible^ tn'C 

fr*gy*zcygfrke scattering by large aeroso!j^.rfj.cJ.§sjn the air and wafer droplets 

in a cloud layer may be significant ./-Swc-s the present UCLA GCM covers/ 

^x/ y th e troposphere, the a b so rptl on_b y ozo n e be neglected >/ Bjrtbermor e r>-H : or 

Aimpli c Ity : ^any s catiec. i ng-o Mo ng^Yave-ra di aHorrrs^TsTiffi^d^^bs^eg'tlgf W e • 
IfptSvXtjfr*’ /em$ w-A I .Scatfo-YWJ U^c// Az -*iGj 

By solving the equation of radiative transfer^ subject to the boundary 
conditions that/T the downward long wave radiation flux asSESsxat the top 
. .. o"PfHe~ atmosphere £§£ the upward long wave radiation flux at 

the earth’s surface Is tisctfcr^Lv^^^lack body radiation at the surface temper- 
ature, the downward and upward radiation fluxes R* and R| at a reference ievei 


z^can be expressed by: 
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R| - / IT B v CT 2 ) dv - J TT B v a ro ) T f - U z ) ] dv 

o o 

03 d B CO 

+ 1 M ~Yr~ T f c \ (u " u 2 ^ dTdv » 

O 7_ 


(H.l) 
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I 
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?A 


A 


*1 •/ 


« » t d B (T). 

R|=X’ r B v g 2 )dv + JJ TT 


o T. 


dT 


VV U z~ u)]dTdv 




(»• 2) 


where u(^) Is the amoun^y n the vertical air column from the earth 1 

surface tc the level ^ , T(j ) -is the temperature t tt B^(T) Is the flux 
of black body radiation of frequency v at temperature T, Is the generalized 
ab sorption c oefficient defined by Elsasser (1942) and Is the transmission 



/ 




I!~2 


s 


function of a slab at frequency v. The subscripts g, z and » denote the earth's 
surface^ the reference level and the top of the atmosphere, respectively- 


Since the width of an absorption line Increases by collision damping, 
the absorption coefficient l is a function of pressure p. This effect can be 
approximately expressed by (p/{^ ) a where p>^ is a standard pressure and a is 
a pressure scaling factor. Since the absorption coefficient Is always multiplied 
by the amount of the absorbing medium u, we can conveniently use the absorption^ 


coefficient at the standard pressure and apply the pressure correction to u. Th 
water vapor amount modified in this manner, the effective wafer vapor amount ^ 


» s ,9 Iven b y 




) J4 


u * := "g^ 9q (V^7^U 
' p( v °° 




/ V 





' (He 3) 


A 


where q is the 'mixing, ra ti o of water vapor d 
used in lieu of u. • • - 


if In the following d 


ng discussion, u is 




Z( u \o i % , r ) s jisy l * ~ar T A, u) dv 


We define two types of weighted mean transmission function as follows: 

(II. 4) 


* 

*2 


dBCT) 


v sr° 
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y , 4 ■ -/ s 


f i 


tbJj f uB v CT)T f ( V )dv ; 


(». 5 ) 


A 


wher 


re 


7tB(T) = JirB (T)dv = oT 


(H.6) 


A 


Relation (11.6) is the Stefan-Boltzmann lav/ and qj the Stefan -Boltzmann constant 
is 0.817 x 10~ 10 ly min" 1 °K' 4 (1.17 x 10~ 7 ly day" 1 °K 


r and t are 

functions of temperature because the spectrum distribution of black body radiation 

1 Iin ^n-I- r. g . — — *-» ■ * " r ~~ ~ 

gradually changes with femperafureA' Yamamoto (1952) calculated the mean trcms~ 
mission function t(u , T) for water vapor and found that: 1) the dependence of t 


zt 


3 



II-3 




on temperature is weak tn the range 210 K to 320 K and 2) below 210 K, 
t decreases very rapidly with decreasing temperature., Therefore, following 
Manabe and Mo Her (1961), if we select some critical temperature T c near 
220°K, t(u , T) may be approximated for T >T^ either by t(u , T ) wifh^ 

T — 26Q°K,ss5* or by a mean transmission function averaged over the^ range , 
210°K to 320°K» 'tfa* r ) ^c(u) Z A 

the calculation procedure for long wave radiation is considerably 
. simplified by use of the mean transmission functions. The second term of the 
right hand side of (11 *1) is rewritten as follows: - •■■■; - 


CO 


i W u « _u z )} dv 

O . •: ; '• • 

>/- B y<V T fMV U .» dv 

co T a> d B \ 

a/;’ '*'* ■ ^ “ ' + U * df Wv u z )ldTdv ; 


A 


A 


V, c-r rv ! : 


V — , * * T \ , A 5 / * * T\ , T 

'= ^3 t(u oTV T c ) + f t(u -“ V T) 11 dT dT > . 

T c.. 

= <*T* r(u*-u* f T c ) + (aT^ - aT^ ) t(u*- u*, T ) . 

Substituting this relation into (II* 1) we obtain 

v ' ; '" ' • R|=ttB z -irB c t(u*-u*, T c ) - (it B^- tt B c ) t(u* - u * , T ) 

irB * : ‘ .* "■ ‘ 

+ ^ t(u - u z , T ) d(irB) , / 


(H.7) 


A 


i: v-':* v 

where 


tiB. h aT: 


/ 


4 
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Similarly, the upward flux at level z Is 


ttB 


* — , 




Rt — it B + f t(u* ~ u 7 ', T) c1(itB) 
z z . J z 

wB z 


The net upward flux at level z becomes 


R s Rt- R4/ = A + B + C . * 
z z z 


(11.3 


(If. 9) 


A 


/i 


where 

< 'i'vr ' 


: A ~ it B t(u - u , T ) . 

C 05 z c / 

^ , lr . 

B = (ttB - ttB ) t(u - u , T ) } . 

x < x > c c° Z / 


ttB 


ttB 


C = J t(u*- u*/ T ) d(uB) + J t(u/ o 7T) d(irB)^ 


g 


* * 


nB w 

ttB 


ttB. 


z 


f 9 T( |0*-0*|,T)d(itB)' 

v ^ 


. A 

A 

A 


(IS. 10)/ 


The relations above are schematically Illustrated by Yamamoto's Radiation Chart 
shown In FIgure(IKK) — //“/# 

The heating rate due to long wave radiation is given by 

_ ... ;• ~ ^ ‘or: • 

: • ' - . ■ ■ - ’ -• • ; •• .r §!.= . 

9t . c 9p ' 

■ . ‘ . p 

v/here g is the acceleration of gravity and c is the specific heat of air at 

P - 

constant pressure * : 


(11.11) i/' 


A 










II -7 



However, t exponentially decreases as ]u - u. | increases and the change 
■ of ju*-u?| with ® Is also roughly exponential . As shown in Figure 11-1 , 

% ' T in the Wighb^+Tood of level i usually decreases more rapidly withyirB than 

a linear relationship. Therefore a careful formulation is required to compute 

the contribution to reference level i from its adjacent layers, C._^. and C. + ^. . 


8 

preceding page blank-: hot. filmed 



To parameterize this contribution, a bulk transmission function t is defined as 


follows; 


' C. .■= T ‘t(u*- u*, T )d(nB) h(uB - irB )t ’ ' 

S n 1 1 1 2 *• 2 /‘ > 

~ -Tf D. •' ' • ” " *• 

1“ 2 


7tB. + 3 

C i+i i = - : J” ■' T(u *“ u * T)d(tfB) a (nB. +a -irB. )t. + ^. 

' -rB. 

i 


Since t - 1 for f u - u. I = 0, t. . (or T. . .) must have some value between 

1 i 1 1“3 ,1 l'T‘2 , t 

] and t. . (or t. , Thus we define a linear interpolation factor m by the 

1-2,1 _ 1 + 2, 1 n .... • 

following relations: s - • 

//'• ;i i _ a/l =;(Lt.«h- 3 ;i V 2 / ( 1 +r V 3 ,i L o (1U8: 

y ■ -V»-,r <f + V 3 ,i W /( 1 + V a ,i> ' ■ (IU9! 

The values of m must be determined by the physical parameters of the 
adjacent layer under consideration . We assume that the vertical distributions 
of the water vapor mixing ratio q and the temperature Tlcan be represented by 

' k. A 

b ! - . q =q i (p/p i ) * ) (U.20 

- ; ’ T=T+r(p r -R) * ‘ (n.2i 

. . _ * » # o 

This, m may be a function of pressure, mixing ratio and temperature at the 
re f erence level i, the lapse rate parameters k. and Xj and the pressure thickness 
Ap of the layer (see Figure II-3). Thus, . 

r , .^(p., q./ T., k,, T., Ap) : - (11.22 


".where m. = m. , ; » and m. = m. . 


i i+3,i 


l l“2,I 


9 



11-9 


W 3 

i \' nJ| 


~ . For the. G CM atmosphere^ in. which the above parameters vary over 

a v/idec^range4han thaf observed in thereal atmosphere /the ihfegralsj.n . 
(11.16) ancUn.,17) are. computed numerically by dividing t he GCM layer under^ 
consideration into thin sublayers of 10 mb t h tc hr^ss«J~ln this way, the values 


P - 


V £" r ", 



? k » ; . 




aP rn u2 g=m 


. - - o cAert ■ftrs 

r>a.j i!~ s ; , ■ , , ■ H8gi*e*8- 

po/* rt*n^ ^UlCh \ *7^ * l^_ r ' J_ ' - - • 

l^/^sestimated from (II. 16) and (i 1.17). Then, from (11.18) and (11.19), 
the values of mf -obtained as a function of p., q./ T., h., T. and A P* 

The functional forrr£of mf then found by a curve fitting technique. 

'•yY yliXv . ^ 

-teimBRaty, calculations have shov/n thaf the dependence of m. on 

T; j\s negligibly small; it is therefore ignored.- Figure 11-4 shows that m. can 
be expressed approximately as a linear function of Ap when the values of the 
remaining parameters are fixed . Thus, 

m? = at + b? Ap/100 . (1L23) 

i i i ■ ^ 

where at and may be functions of p., q./ k« and I\. Since the dependence 

of a* on k. and XI is weak, it is neglected and the following empirical relations 

I i « 

are obtained: ■ . ■ 
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L + (P.)+ F + (Z.) , 
a i a i 


s 


Ar -,- 

Lq (P j) = "1-66 + 1 .76 log ]Q P., 

(P.) = -0.197 + 0.0002 P. , 

F + (Z.) = 0.30 Z. + 0.28 zi> + 0.04 Z$, 
a i . i i i 

F^(Z.) = 0.0812 Z 4 - 0.045 Zf>+ 0.02334 zf, 

* ' . 

= Min (-0.041 + 0.021 Z., -0.06 )< 0 , 
(|| ;) + = Max (0.01225 + 0.007 Z., 0.0093) >0 . 

-0.09 L" (P.) + F - (Z. - 0.105 L"fP.U . 
a j a * i a v i' ' * 


( 11 . 24 ) 
( 11*2 5 ) 




{ll. 

( 11 — 26 ) 


-0.0? L- (P.) , ♦ F- (Z. - 0.105 L- (P.» * (I)- At. ♦ (II)- AT. . 01.27) 

L" (P.) =Max (61.86 - 22.92 log 10 P., 76.63 - 28.39 log P.) , 
L~(P.)=Min (-42.59+ 15.78 log^P., -60.81 + 22.53 log P.) , 


F (X) = 2.57+ 0.233 X+ 0.18 X3>+ 0.027X3>, 

a 

F~(X)= 1.42 + 0.43' X + 0.16x1+ 0.011 xJ>, 

.At * * 

(§kj " °’ 08 + .(0.371 -0,102 log 1Q P.) (Z. + 2.1) > 0 , 


i i 


(|— ) = Min (-0.0325 - 0.005 Z., -0.0275) < 0 , 




' 0 


12 
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where Z. - log.jQq., Z- 1 - |Zj+2.5|, Ak. k.~3, 

Ar. = r.-10 (°K/1QO mb) and X is a dummy variable . 

Figure 11-4 shows the relation between and Ap for the cases of 

q. = 10, 1,0.1 g/Kg, p. = 600 mb, k. = 3 and I\ = 6°K/l00 mb. The 

plotted points give the values of mf obtained from equations (H£16) - (IK19 ) ^ 

while the straight lines represent determined from (11.23) -(11.27). The 

agreement eppecN¥*to4»e satisfactory. 

A 


& 


The values of a? and b^ for various values of p. and q. are pre~ 




sented in Table 11-1 for the case of k. ” 3/ .32 _ 6 K/100 mb. a? and ^ 
increase while b. decreases as p. Increases. Their dependence on q. is not 
monofonic. a? gradually decrease with decreasing mixing ratio and reach 

l i 

minimum values near q. ^ 0.1 g/Kg. b. exhibit a similar behavior with 

1 1 + 

minimum values near q. = 0.01 g/Kg. Since b. is considerably smaller than 

* * — 
bo , while the difference between a. and a. is usually less ihan 1/ m. Is 

1 + 1 . 1 
larger than m. in most cases. 


~v. 




r-the 




As discussed in Appendix C, if the permissible error irnv^/" .- and 


G.Z7^\is 0.2 °C/day, the accuracy requirement for m. is not too severe and 

r \±± iy * ± 

usually more than 2 0^ error in m. -co^ld^be tolerated. 

If the vertical coordinate of a GCM is pressure or geopotential height, 

p. and Ap are exactly or nearly constant In space and time. For such a GCM, 

I . i . . 

only q., k. and Ij are variable parameters in the expression for m. . Considering 

the range of these parameters in the actual atmosphere, appropriate mean values 

for k. and 32 , such as their globally averaged climatological values, may be 

used in the empirical expressions for m-. For example, k =“ 3 and T ** 6 K/100 mb 

may be adopted. • 
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TABLE IM * The values of a* and b* for various p. and q.. The values 












IM4 


/'K/ 

hit- 

f- 


As a rough approximation, the effect of the variation of q. may ' 

also be ignored. In that case, the hemispheric mean climatological value 
given by (IV* 16) is recommended as an appropriate constant value of q. 

However, the mixing ratio at a reference level In the lower troposphere varies 
by more than two orders of magnitude over the global domain. Therefore, it 
is desirable to retain the spatial and temporal variation of q for the estimation 
of m?* 

i 

When the cr-coordinate system is employed^ in the UCLA GCM, 
p. on a g-!evel and Ap between a- levels may van/ in space and time. For 
example, over high mountain areas such as the Himalayas and the Antarctic 
continent, p. and Ap are significantly reduced^ which leads to a decrease in 
m? and an increase in t*. If p. is assumed constant on a reference a-leve! , the 
contribution to the flux of long wave radiation at the reference level from its 
adjacent layers Is underestimated over high mountain areas* 

^3.,; The effect of clouds ; .• 

/ ' . ’ . • 

When the sky Is covered by a thick cloud, the top and bottom of the 
cloud treated as perfect black body radiators at their respective temperatures 
(the net flux In the cloud is assumed to be zero). This assumption can safely be 

°^a<£&ee{-for a dense cloud.buhmust be checked for cirriform clouds which contain 

A ' S 

small amounts of wafer vapor and solid wafer (ice crystal concentration^! ess than 

Oo Vg-m ^)* • " 


Friiz and Winston (1962) analyzed the radiation measurements of TIROS II 
and found that the cirrostratus layer did not strongly absor£ythe window region of 
•te long wave radiation, 8 p ~ 12 p* A detailed discussion of the greyness of 
cirriform clouds has been presented by Katayama (1966)* 
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/ ' do, o+ ' - . -***&&& ^ 

ilrrrform clouds^completely absorb long YiQ^ md\aUor\j simplicity vi**- 

as grey body 

radiators with greyness of 0.5. In the computation of the radiation flux, the 
above assumption ^-be introduced roughly by modifying the amount of cirriform 
clouds by the weighting factor of 0.5. 

It has been confirmed by laboratory experiments that there is a critical 
temperature , about -40°C , at which homogeneous nucleatlon transforms all 
liquid water to ice crystals. Clodman (1957) found from aircraft observations 
that the base of cirrus clouds is situated at the level where the air temperature 
is about ~40°C . We therefore assume that any non-convective cloud formed in 
a layer with temperature less than -40°C Is cirriform and that the greyness of 
the cloud Is 0.5 


1 — — ’ ,.. T .,-. 1rn -,.. ^ 

. // /— 

' ; 4; Empirical transmission function equations ^ 

y / 1 , 1 , . .. . rr __ 

a - / — * — ^cjnr^Tf 7 ! . 

v " p ' f AG^pd 4 ^g->ro^amamoto (1952), the total transmission function of a 

/ mixture of water vapor and carbon dioxide is assumed to be represented by the 
J product of their respective transmission functions/ that is 


-^Xo^co/ 10 “ t h 2 o (u h o^ t)t co (u co J) * (!K28) 

0 Kl^ ’ 1 v 2 • . y £ * * *• A 

■Mr' (> ^ uVmg the tuple presented by Yamamoto (1952), we deduce the \ • - 

Y'p 'Jr * t . i 

pM hj^/ following empirical equations for the transmission functions of water vapor 


J and carbon dioxide . 


C /! ' 


■ , ,; T . .V, /V// 


(i) Water vapor - — - — 

— -j Two empirical equations for the transmission function of water vapor at 

/a temperature of 260°K are presented below; the first has a form similar to 

\ Calendar's expression (1941) for carbon dioxide while the second is expressed 

\ 


/a , t d> ' / 


, \ 


/ / J r \S jV-* ‘tl* ’ 


/,,* T-)> A/*' 7 


Si 


by a quadratic polynomial 


(A) T o^ u * 260 ° K ' ) 1 *0.416 

h 2 U 1+1.750 u 


(II .29) 


r'-T" * 

.(B) T (u* 260°K) = 0.384 - 0.2785 Z + 0.0345 Z& (u zl) ; 

; • • = 0.384 - 0.259 Z - 0.028 Z®, (10 -4 =5 u*<l), (11.30) 

= 1/(1+288 u*), (u*<10 -4 ) / 

ie . it "k . 

where Z = Iog,g u an ® u ” U H O 0 

. 2 A 

' Table 11-2 shows the accuracy of the two expressions above . In the real 
atmosphere, u* is usually less than 6 g cm . Therefore the maximum error 

of expression A is about 0.025 while expression B has a negligibly small error 

; . ; ^ —2 "2 
for u between 10 g cm and 0.0001 gem . 

1 1 -* 2 . * .... e * 

TABLE Accuracy of.empirica! equations for the transmission function 

O 

of water vapor at a temperature of 260 K. The errors are 

S "3 ’ 

shown in unit. of 10 . 




0 

6 

3 

1 : 
6.3 
o;i r 
0.03 
0.01 
0.00 
0.00 


0.140 

0.183 

0.259 

0.384 

0.510 

0.612 

0.713 

0.790 

0.910 

0.972 


0.180 

0.212 

0.266 

0.364 

0.485 

-0.598 

0.711 

0.795 

0.910 

0.963 


0.140 

0.188 

0.259 

0.334 

0.512 

0.615 

0.714 

0.790 

0.909 

0.972 


— ■> 
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TABLE Accuracy of the empirical equation^for the transmission function 

of water vapor averaged over the temperature range between 

, - 220^^}300°K. Jl\3l 


* , . . 

u H 2 ° 

(g cm 2 ) 

r • . 

* - 

A t * \ 

" rH 2° UH 2 0 




Obs. - 

C 

o-c 

• D 

O-D 

E 

O-E 

10 - 

0.134 

: 0.172 

-28 

0.134 

0 

0.134/ 

0 

6 " 

0.181 

0.204 

-23 

0.185 

-4 

0.181 

0 

3 

0.250 . 

. 0.268 

-18 

0.248 

2 

0.250 

0 

i ' ■ ■: 

0.373 

0.364 

-9 

0.373 

0 

0.373 

0 

• 0.3 

0.499 

0.489 - 

10 

0.498 

1 

0.501 

-2 

. .0.1 

0.603 

0.597 

6 

0.601 

2 

0.605 / 

-2 

0.03 

0.704 

0.715 

-11 

0.704 ' 

0 

0.705 

-1 

••• o.oi ■ ■••• 

0.783 

’ '0.800 

-17 

0.784 

-1 

0.782 

1 

0.001 

0.905 ‘ 

0.914 

-9 ’ 

0.905 

0 

. 0.904 

1 

. -0.0001 
t 

0.971 

— — 

■ .0.9 66 . 

5 

0.969 

2 

0.971/ 

0 


As previously mentioned, the transmission function of a slab of water vapor 

is- nearly constant with respect to temperature; it varies by about 0.02 from 220 K 
to 300° K. it may therefore be more suitable to use an averaged value of the 
transmission function over this temperature range rather than the value at 260 K. 

A 

We define q as * • - 

' : •• 2 \ rs ) * C:.2uf 

4o<"*> ■ KT„ 2 ifa*. 30O°K) +:t^o ( u * 2«o“K) t 22 ‘ fK)] 


(11.31} 


The ^se^ ej^alues of q are shown in the second column of Table Il“3. 
From the observed values, 2 three different empirical expressions are deduced 
as follows s 




IM8 


y 


( c ) Q (u*j - 


1 


1+1.746 u 


*0.423 


(Dj Q (u*) - F(1 .681 ,0.535), (u"i 0.724)^ 

— - — = F(1.600, 0.3815), (0.00186 ^u*< 0.724) 


J 


* 


= F (3.853, 0.522), (u < 0.00186) i 


where F(a ,b) - l/(l+au _^)« j \'__i 

(E) -A n (u*) = 0.373 - 0.274 z + 0.035 Zg) (u"* 1). 

. :=. 0.373 - 0.2595 Z - 0.0275 22) (10~ 4 s u*< 1} 


i/jf+ 298 . 7 »*) t (u*< io" 4 ) , V 

. Expressions (C) and (E) are similar to (A) .and (B), respectively; expression (D) 
is an extension of (C). The accuracy of these three expressions Is shown In 
Table !1~3. Expression (C) exhibits the largest error, but its maximum error is 
~less~tban 0.025 for u ^6 gem «- „ 


J 


(11,32) f 

' ‘ ' 'A 

(H.33) A 

* 

A 

A 

> yj 

(11.34) • A 


'~V* 


It is also necessary to evaluate *t(u , T), the mean transmission function 
'defined by (11.5), As suggested in section II. 1, 220°K is adopted as the critical 
temperature T . The empirical expressions for t are similar to those for t, namely 


W T H 2 0 (u ' 220 1+3 0 u *0.403 


J: 


ti 

\ 


, *2,0-1 (|[.35) 


fa. 


* 


A 


(B 1 ) ? H Q (u , 220°K) = 0.254 - 0.1985 Z + 0.0205 Z2) (( 0.1 s u < 10^ 

= 0.216 - 0.2827 Z - 0.0253 22), (lo” 4 £u*<0.1). 




(11.36) 


■ . = l/(l+2.56 u* 0 ' 39 ), • (o*< 10" 4 ) V 


y l 


L\*l\ 


A 

A 

A 

A 


'(C) Q (u*, 220° K) = 0.254 - 0.2224 Z + 0.0444 zl> (fa £ u*<K)/ y A 


1 


= 0.254 - 0.226 Z - 0.007 Z3i (0.03 *u<l) (11.37) » A 

= 0.194 - 0.297 Z - 0.028 23) 4'0$(v*< 0.03), % X 
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A 

A 


y 


The accuracy of these empirical expressions is shov/n in Table 11-4. i 

Expressions (A) and (A 1 ) were used for q and q/ respectively, . 

in the 2~Iayer moist and early-stage 3-layer UCLA general circulation models. 

; \JC LA 

In the cun*ent 3-layer^GCM, expressions (E) and (B 1 ) have been adopted for 
Tj_i q and respectively. ... . - A ... ‘ 

— n-y, - 

TABLE Accuracy of the empirical equations for q at the temperature 

(i •/ , cf220°K. The errors are shown in unit of ‘ 10 


* 

U H 2 ° 



T h 2 ° 

(u*^, 220° K) 

• 


-(g cm 2 ) 

- -Obs. 

- : ' A 1 ■ 

O-A' 

... C‘ : 

" OC 

B‘ 

O-B' 

10 . , 

• 0.076, . 

0.115 

. -39 

0.076 

0 . 

0.076 

0 

3 * 

“0.158 

0.176 

-18 

0.158 

0 

0.164 

-6 

1 — 

-0.254 - 

0.250 

4 

0.254 

0 

0.254 

0 

0.3 

0.368 

0.353 

15 

0.370 

-2 

0.366 

2 

0.1 

. .0.473, . 

0.460 

. 13 

0.473 

.... 0 . 

.0.473 

0 

0.03 

0.584 

0.582 

“ 2 

0.582 

2 

0.585 

-1 

0.01 

0.676 - 

0.685 

-9 

0.676 

0 

0.678 

-2 

.0.001 

0.833 . 

0.848 • 

-15 

0.833 

0 

0.832 

1 

0.0001 

0.934 

0.935 

1 

0.934 

0 

0.934 

0 

0.00001 

0.972 

0.973 

1 

0.979 

-7 

0.935 

-13 


- (i!) Carbon dioxide " — ; " '' ‘ : ' ' 

The transmission function of carbon dioxide at norma! temperature and 
pressure (NTP) Is shown in Table 1 1 -5. These observations can be expressed by 

r C0 2 {ii C0^ = 0.930 - 0.066 \og }0 u* Q ^ - <•}■?'- 

•A .• A,-. - .for u^q £ 6.3 cm - NTP^ (II *38) 
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A 


>v- 


where u ^q Is the effective amount of carbon dioxide at NTP within a layer 
of thickness^ | p. — p, j (see section IV .2, equation (IV.19)). The agreement 
between (11.38) and the observations is good^as shown In Table H-5* Introducing 
(IV. 1 9) Into (II. 38) yields 


A 


t co (M> = 0.791-0.066log 10 [ 


1 L 

r\ 2 


oo 




(11,39) 


A 


The temperature dependence of the transmission function of carbon dioxide is 
neglected in the present formulation. 

-• 1 1-57 

TABLE IfcF. The transmission function of carbon dioxide. Cal^shows the 
. ol>h\AU 

value^si^a'ted-oy. v*» — — 



J^ 0 ^(cm/NTP) 

1 

3 

10 

30 

100 

150 


Obs. 

0.930 

0.900 

0.865 

0.834 

0.796 

0.734 

t co 2 

Calc, 

0.930 

0.899 

0.864 

0.833 

0.798 

0.786 


If pressure is the vertical coordinate of a GCM, ('*0 does not 

change in space or time. In the a-coordinate system^ however, wee-4>e-* 

^ ca ^^retsure at a constant a-levei varies with the surface pressure^, t cq (i,j) 

also varies with p . Thus, with the exception of the layer above the a ^ 0 
S , ->nu£ir ^ 


A 


T7~ $**** 

level, T_ n (i,j) ske**S3 be. calculated at each' grid point and<$? each time step. 

From an examination of the dependence of (i, j) on p^ the following simplified 

empirical relation is adopted: ■ 

(11.40) 

where \ p = l+0.09(1-p $ / Po<) )^ V (H.4I) 


t C o (Ui;p s ) P T cod' ,i;p <V 






A 

A 


A 

A 
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and T c02 (i,j; j|) is t co (i,i) when the surface pressure is p g and p oQ - 1000 mt. 
Figure 1 1-5 shows the values of t co jfor various layers of the 3-layer GCM. 

, 

Finally, the total transmission function for the layer between levels i and 
j Is given by . 

V s u co 2 ( ' , ’ ) ' T) ; 

(II. 42) 


where 


l 2 w "2 
^S } T C0 0 (I ' i; O W M) J 


H~0 


(M) = t h 0 ( I^V) ' 


5. Finite difference schemes 

— :_v j n t his section, the finite difference scheme for the long wave radiation 

calculation is summarized. The indexing scheme is shown In Fig. Il~2. The 
following assumptions are made regarding cloud layers: 

... V - . .... . • 

(i) a cloud layer is located about a level with an odd index 
• (ii) the cloud thickness is equal to the thickness of the layer in 

' ’ * F which It is situated 

(iii) the distribution of cloud elements Is horizontally 
quasi -homogeneous 

' 1 (iv) there are no clouds in the layer above level 0 

(v) the cloudiness and greyness of a cloud in layer j (between 
; • r / levels j-1 and j+1) are CL. and a!, respectively; the 

modified cloudiness CL*, is defined by CL . - a. CL.« 

Combining (11.9) and (11.13), the net flux of long wave radiation at 
reference level i in a clear sky is 






In a cloudy sky, (11.43) is modified as ■ 

l-i j- 1 

R. = [irB t . +(ttB - tt B )t . +C . +C .] tt" (1 -CL ‘ ) +C 3 ; .n (1-CL')+.. 

J *- c «>» . » C. «l -1,1 1 / 1 j-- ^ j * j=3 I 

, ‘ % A-Sn’ c-awo^ 


1-3 

+ . » . + C . '(• it. 1 .- (1-CLJ) 

I" 3/ 1 j— i+i I 

+ [C_ , . + (wB - ttB )t ] VO-CLj) 5 
i-i/i g 1 I* j=i+i l y 

where tt! denotes multiplication over j odd. t.. is expressed by (11.42) 


(11.44) 


Substituting (11.14) and (11.16) -(11.19) info the above relation yields 

R. = YcrT 4 t . +AB t . + AB 
i \ c c, ro 03 1 «v 


T . + T . 
«l Ot 


T .+ T • l-J 

I o l \ * 


+ AB °V -) »’ (1-CU)-.. 

0,2 2 y j =1 J 


✓ 


/ 


+ AB 


T s ! +T 4i *“i 


2 ,* 2 


TT* (1CL!) +. . . 
j-3 I 




1+m. t„ 

1 1-2,1 


i-a,i 

l+rrC 

i 


+ 

1 + m. t . 

1 1 + 2,1 

i,i+s 

~ l + m+ 

. 1 


T. .+ T. 

1 + 3,1 1 

i+2,l+4 2 


l-l 


i+i 
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+ AB i-^ . i 1-3 v> (i-CL!) 

1—4, 1 — 3 4 j—i+i 1 




•+ AB t 

1/9 I 


I-i 


/ty/ i=i+i I s 



where AB - crT 4 ~ crT* . 

A,m rn A - 


As an example, the finite difference scheme for the current UCIA 
3-layer GCM Is shown In Appendix A. 7 
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III. SOUR RADIATION 

There are three principal factors that contribute to the depletion of 
solar radiation in the cloudless troposphere: 


(i) absorption by water vapor 
(ii) Rayleigh scattering by air molecules 
(iit) absorption and Mie scattering by aerosol particles 


According to Katayama's preliminary result (1966), absorption by 
aerosol particles is comparable to absorption by a cloud layer In the normal 


mean cloud state and Is of the same order as the absr^ption by water vapor, in 

i high latitudes at large sun zenith angles. Thus the depletion by aerosol 

absorption should be included In the atmospheric radiation process. At present 

however about the optical properties of aerosol particles 
, ^ ' 


A 



, ; V 




A 


V, Basic quantities 

The effective absorption bands of water vapor for the solar spectrum 
are in the wavelength range X>0.9 p . There are two absorption bands in the • 
region X < 0.9 p, namely \= 0.75 p and 0,85 p, hov/ever their contribution 
to the total wafer vapor absorption is negligibly small* 

On the other hand, since the amount of Rayleigh scattering is Inversely 
proportlcnaf to the fourth power of the wavelength X/ the scattering by air 
molecules rapidly decreases with increasing wavelength. With an error of a few 
per-cent, scattering in the v/avelength range X> 0.9 p can be neglected, 

(tucjiqyp) i 

Based upon these considerations, Joseph divided the solar 

radiation into two the scattered and the absorbed parts. This idea 

A . . - 



26 


was also presented by FeigePson (1964), using a critical wavelength of 0.75 p 
Instead of 0,9 p. Thus,for X<0.9p, solar radiation Is assumed to be Rayleigh 
scattered but not absorbed by water vapory while for X> 0.9 p, the converse is 
assumed. The solar radiation incident at the top of the atmosphere is then 
separated as follows: - 

The scattered part: = 0.651 S Q co$£ (Ilf . I) 

The absorbed part: $£ = 0.349 $ Q cos£ (111.2) 

where is the solar constant and £ is the zenith angle of the sun. 

(i) Absorptivity of water vapor 

Mugge and Moller (1932) found the following empirical relation for 
the absorption by water vapor of the mean solar constant: 

. .. a - 0. 172 (u sec£) ly mm * 

where u is the effective wafer vapor amount. ' . 


Yamamoto and Onish! (1948) and MacDonald (1960) also constructed 
absorption curves. Their values are smaller than those of Mugge and Moller 
by about 15 percent, Manabe and Moller (1961) recalculated the absorption 
and obtained a 10 percent larger value than Mugge and Moller. 


Since absorption by aerosol particles Is neglected In the present treat- 
ment, Manabe and Moller * 5 larger value of the absorption by water vapor^ 

o‘ s 0J8? (u sec^)^*^ ly min' ^ (III. 3) 

is employed. Since this absorption refers to the mean solar constant S q , the 
absorptivity for the absorbed part is . 


a 


•_ ^ / * „\0.303 

“0.271 (u sec£) 

0.349S ^ J 


( 111 * 4 ) 


where S = 2 ly min 



0 ! 




,o'V 


m-3 


,0.303 


Defining theftunctlon A(x) by 

A(x) = 0.271 (x) v ; 

the absorption by wafer vapor of the direct solar radiation is 


(III. 5) 


a 1 - S° A(u sec Q 


) 


v/here 



= 0.349 S q cos £A(u^ sec z) 





(!». 6 ) 
(in. 7) 


and r- is the sun’s distance from the earth and r £ Is the mean value of r E , 


that Is one astronomical unit. The zenith angle £ Is given by 

cos£“sinc p sin 6 + cos<p cos 5 cosh^ (U1.8) 

where <p is the latitude, 6 Is the solar declination and h Is the hour angle 
of the sun. . ' * . = . 


In a numerical simulation of the seasonal variation of the general 

fc ->i o w *>i * 

circulation, the seasonal variation of 6 and r e must be f>red+eted^rd-abw4a4-ed'. 
As shown in Appendix 0,^6 and r E can* be -estimated by a perturbation solution 
of Kepler’s second law. . — • -*' u v ' ** 0 ^ ^ 

■ t \J&- 

(M) The albedo of the/bfm^sph^e^ ^ 

-=■ ' The albedo of thc^ctacff-sky due to Rayleigh sea tier I ng^ was estimated 

by Coulson (1959) and is presented in Table 11-3 of the paper by Joseph (1966). 

Josepl^t&d-k^xpfe&s the albedo by a least-squares polynomial in cos £^of the 

form • .*• 

a o (?) = S C. cos' £ » 
i 

In the following, the empirical expression 


.1 


P. 


«j 


a o fe, p s ) = 0.085 - 0.24§)log 10 M cos 

oo' 

. will be employee^ where Is 1000 mb and p $ is the surface pressure, 

r a Table IIM gives a comparison of (111.9) for p "p with the values of 

Coulson, The error of (III. 9) Is less than 10 per cent* 


(Hi. 9) 


A 


A 


A 


A 


y\ 


A 


A 
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(Ht) The albedo and absorptivity of clouds * 

The albedo and absorpfi vlty of a cloud are functions of the cloud 

thickness/ cloud height, liquid and/or solid wafer content, wafer vapor content 

and zenith angle of the' sun. Since observations are sparse and theoretical 

studies are Insufficient, it is not possible to introduce the above effects explicitly 

into the model* Therefore, following Rogers (19 67 ), the albedo and absorptivity 

of clouds are prescribed as shown in Table !!!~2 (here subscript i denotes cloud 

type)* The definitions of u ' 1 and u will be given later. 

cr ct 
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?» Interaction between cloud layers and solar radiation 

A o ftfic Softer ; 

Por the scattered molecular scattering in 

clouds is negligible compared to scattering by water droplets. Absorption by - 

water droplets and water vapor in clouds can be neglected. 

For the absorbed part scattering by wafer 

droplets and absorption by water vapor and droplets in clouds must both be ■ 
considered. Thus the behavior of the absorbed part in a clqud layer is very 
complicated. As yet, an analytical solution suitable for the present purpose 
has not been obtained. Therefore, the absorption and scattering of the absorbed 
part in a cloud layer are treated separately under the assumptions presented 

in subsection (Hi)* •' - 

ftn 


&& , 

(i) The albedo of a cloudy atmosphere^scattered parfj. 





Consider a single cloud layer with albedo R‘ for ihe scattered part. 
-Since absorption by the cloud is assumed negligible for the scattered part/ the 
transmissivity of the cloud is (1-R‘). The transmissivity of a cloudless 
atmospheric_co!umn due to Rayleigh scattering is (1-qJ . Thus the fronsmissivity 
^for the scatteredjra i^oT^ re with a single cloud layer/^an be roughly 

estimated by (1 — R ' )(1 * Then the^ajbe^ 

atmosphere with a single cloud layers 

; ; ► a =1 -(H’)(W ) ; 

- C O' 


= R* + a - R‘a 


(HI. 10) 

'•f 

o o 

For a partly cloudy sky v/!th cloudiness CL, is modified as follows: 


a -(l<:L)a o + CL-(R‘+a o ~R ! a o ) } 


a + CR’(l-a ) . 

o o / 


CIII.11> 

where CR' = CL-R'. (ill. II) is identical to (111.10) with R' replaced byCR'. 


A 

o 

A 


,/A 
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(H) Tne albedo of multi-layer clouds ' ' ‘ ' ^ 

(a) 2 cloud layers . . h 

Let the albedos of the two cloud layers be R l and . Multiple 

reflection between the two cloud layers must be considered „ As shown 

schematically in Fig- 111“!/ the total transmissivity is the sum of 
T .I T II T III , 

*13* *is^ 'i 2 /“•••/ where , 

; fi 3 = o <)(i<) , 

N T S- = o-r;)(i-r s )r;r; , 'y 

J T ” 1 = (l-R;)(l-R,)R»Rg . 




0-R, 


(i-r;xi-r2) 


12 - i 12 . 


F' tcj, //A A F 77 c <£r$\a rt ySA AA. 

-tohJ jyjo 




HI-7 


Then 


4 


./ 


- 


13 


T J + t ii + T m + 

'ia *13 *12 


R 


(I-R x )(l-R a )ci +M 3 + fl^R ,) 3 + 
0-r;)(1-r' 3 )/(1-r'r') J 
0<)(1<) 


13 


5 -t 13 = i 




(111,12) 

(ill - 13) 


/v 


/\ 

'a 


( (b) 3 cloud layers 

V:;/ y ; 7 . 

' - Let the albedos of the three cloud layers be R x , R 3 and R i3 . The 
total albedo of the 3 cloud layers can be found by considering two cloud 
layers with albedos R 12 and R 3 • Then from (III « 13), 

^123 “^123 J ' • '• '• • • 


= 1 - 


= 1 


0-O(i-R3> 


13 


; 


(i-r;)(i-r 2 )0-r 3 ) 


(111.14) 


A 


A 


/l 


.• \ 
t \\ 

l fj 


•/ V 


HR^+R-jRa+fyg + 2 R 1 R 2 R 3 

(c) Multiple cloud layers 

, ^ 

Let the albedos of the n cloud layers be R l7 R 3 , 

Repeating the procedure by which (111.14) was obtained, the following 
general expression for the total albedo of the n cloud layers Is obtained: 

s / n / . - 

where E(R.,... / R ) - tt (1-R.) . 

1 n ' i ) ». 




A 


D(R^,. oe/ R ) = 1 - Z R^+2 Z 

i< i * i<j <k 1 


-3 Z 


\<\<k<& 
n~l 


R.R.R. R + , * • + 
i j k i 


(n-l)R 1 R 3 « « * R^ j 


A 
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and, for examole, £ R.R. means the summation of all R.R. for which i < i 


/ £ n i <! 

* r -*' 

and Jfhe number of terms Is equal 

^(nj^elqul valent cloud water vapor amount and assumptions for the 
absorbed part ^ rW / Vf ' ,r»~ 


' >1 



'^4— . '/j • ,j t 

The following assumptions are made for the absorbed part f *■’ . 

(a) The albedo R. and absorptivity A. are prescribed as shown in 
• Table HI-2. 

” / (b) The albedo of a cloud for the absorbed part results from back \ 

I ’ scattering by cloud water droplets. This reflected insolation __ ^ 

■ \ may, before it exits the cloud top, be subject to considerable 

^-/absorption; however this absorption is neglected. . 

, ijf ( c ) The absorption of the reflected radiation after it leaves the 

•. cloud Is neglected., • ; ■ 

| ( j (d) Solar radiation in a cloud and radiation transmitted by a 

cloud are diffuse radiation. The effective optical thickness 
’ for diffuse radiation is assumed to be 1 .66 times the length 
of a vertical column, that Is sec ? = 1 .66 (? = 53 ). 

hjj (e) The cloud layer is regarded as a fictitious water vapor layer 

• . ' • vvhich absorbs the same amount of insolation as that absorbed 
, by the cloud. The water vapor amount of such a fictitious 
• Invnr Is denoted by u^., the equivalent cloud v/ater vapor 

... - amount . 

Let the effective water vapor amount between the top of the atmosphere and 

the cloud top be u* 1 . Since the absorbed part of the solar radiation incident 

on the cloud top Is S a [ 1 - A(u* t ‘ sec?)], the total absorption by the cloud is 

A $° rl~A(u*‘secr) 1? By assumptions, (b) nf&(k), this total absorption is 
i o L ct b 

' "■ ' ' j . 
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£>/ 


equal to the absorption^? a fictitious layer of water vapor having an equivalent 

^ “k 1 It J 

cloud water vapor amount u . . By assumption (d), the following reajf ions^S^^ - 
is obtained: . ^ . • 


* 


A. S [1-A(u 1 sec £ )] 

I O CT 


(1~R.) $° [A(u*' sec T + 1 .66 u . ) 
i o L ct s ci 

“ A ( u * t ' sec ^)] 




By (Ulo5), the above relation yields the solution for u 

* r/A.tl-AWl+O^JAM, 3 - 3 n/ 
u ci LI '0.271 (I-Rjj ; Y 


1 .66 


i ' ' *, 

where x = u seer 
ct s 


(111.16) 


/t 


) (HU 7 ) A 


— irT - 

Fig. UL 2 shows the relation between u . and x for the four com- 

C! * 

blnations of R. and A. in Table Ul-2. The values of u . shown in fable IH-2 
1 * * C1 * 

correspond to the possible range of u ‘ with sec£ = 1 . The values of u . 

* ' y ^ f q ( / /"* *2^ C T ^ * 

('within parentheses represent the mean values adopted in the current model. 

• ' • A *V\ 





The above procedure for the estimation of u . is somewhat sophisticated. 


lues o 


Xu* in Table lifr-2 show'only the normal 1 conditions for the four 


types, 
a cloud t 


y" ; \ \ _\ 

the value of u\ is fixed for each cloud type, the absorptivity 

C \ L *1 'V . \ ^ 

anqe with u ; ah increase or u ^ * 

\ r , CT ^ \ . c t 



^correspo 


onds to a 


In the current model, a cloud with a temperature less than -40 C is 

assumed to be a cirrus type cloud. The normal value of u . for cirrus is very 

C I 


-2 


small, say 0.03 g cm . Thus, when a middle or low cloud is assumed to be 

cirrus, u . may be smaller than the water vapor amount of the environmental 
a • ^ 

layer. To avoid such an unlikely condition, the fixed value of u^. for cirrus 
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is modified by adding to if the water vapor amount of the clear area in the 
corresponding layer.- ' . . •- . 


3. The finite difference scheme fora clear sky .• . • 

The indexing scheme Is shown in Fig. Il“2. In the following/ super - 
script o denotes conditions for a clear sky. 


The downward flux at level l of the absorbed part of the solar 
radiation is 


S°. =S o n-A((u*-uf) S e °£)3 

a i o * 

011-18) 

A 

where Index T denotes the top of the atmosphere. The absorption of solar 
radiation in the layer between j~l and ]+l , AS? , is therefore 

• . 


A C° - C° C° 

AS. =S . -S o 

1 api a j+i 

(111.19) 

* 

A 

. The solar radiation reaching the earth's surface is: 



The absorbed part: $° = S P “A(u T sec £) ] - 

(111.20) 


The scattered parr: S^. = S*(l ” a 0 )/(l“ a 0 a s ) j 

(111.21) 

• A 


where a is the surface albedo and the denominator for the scattered part 

s 

represents the correction factor due to multiple reflection between the 
atmosphere and the earth’s surface. The total solar radiation absorbed 


by the earth's surface is then 


S° - (l-a)[S°+S° T ] # 
I x s L al si 





HM2 


s ' - • 

4. The finite difference scheme for cloudy sky . 1 

(i) One cloud layer 

When the sky is covered by only a single cloud layer, the complex 
behavior of the solar radiation can be calculated by use of the concept of the 
equivalent cloud wafer vapor amount. In the following, superscript c denotes 
conditions for a cloudy sky. v 

From the assumptions discussed in sections 2(i) and (Hi), the down- 
ward flux of the absorbed part of the solar radiation is (see Fig. 1 1 1-3): 

(a) at a level above the cloud 

V. S a U = S o a-A((u T *-u*) S eci;)] } ' ,'(l».23) 


(b) at a level within the cloud 



) 
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\ P \ / 

/ (c) at a level beneath the cloud • 

// S= l =S°(1-R c )[1-A(X+1.66(u* + u* 8 -u*))] ; 

(d) at the cloud top 

s c = s a o - a(x) ] , ;v ■ 

qct o L ) 

where X s (u T * - u* )sec£ . Subscripts U,,M and L denote the levels 

4rn> o-fj Ct i\ . .. « i i 

above, within, and beneath the cloud, respectively. Ap . -is the pressure 

'tHTc!<hesi^of the^Jayer between CT and M. R - and are used in lieu of 

, 

R. and u respectively. , 

The absorption of solar radiation in the layer between j-1 and j+1. 


is then; 


(a) if both levels j-1 and j+1 are above or beneath the cloud top. 


v as? = s c ; ~s c .. j < 

y : ' j a a |+i y - 

: ' > , (b) if level j-1 is above, and level j+1 beneath, the cloud top, 

AS? = S c . -(S°.. +RS' ) . - .. I 

I aj -1 Q |+i J acr .. - i; .\ 

: L{ The last term on the right hand side represents the radiation 

• U ^ d— - of the absorbed part reflected at the cloud top. 

‘ r ■ v • 

The solar radiation reaching the earth’s surface is divided into the 

absorbed part S C „. and the scattered part S C . The absorbed part is given by 
al si 

(111.25) with u* = 0. The scattered part is given by (111.21) with replaced 

by a .where by (HI, 10), 

■ cy s 

a = Hl-R )(l-a ) . 

c c o ' 

i - ' * ’ , 

1 ■ / Tfie total solar radiation absorbed by the earth's surface is then 

S? = (l-a)[S C •! S C T ] , 

, I x sal si y 

' * 
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where 


S UI = S°(l-R c )[l-A(X + !.66(u* +u *))]/(l-R c a s ) ;f 

S sl = S o (,_a c )/(, ' a c a S ) ' 

The denominator of (HI. 30) Is the correction factor due to the multiple 
reflection between the cloud and the earth's surface* 


v- (111*30) 4 

(II I o 31) ^ 


For a partly cloudy condition with cloudiness CL, the absorption :i 

In a reference layer AS. and at the earth's surface S^_ are approximated by 

AS. = (1-CL)AS? + CL’Asf . (Illo32) 

I l I 

\ S x = (1-CL)S° I +CL-S= (111.33) 

K 3V?— <£*/ 

The above single cloud layer mode! in the UCLA 

A 

moist 2-layer general circulation model. detailed description of the 
programming Y? presented by Gates, et al (1 971 ) 0 

(i!) Multiple cloud layers . . 

The assumptions presented in section 11.5 regarding the distribution 
of cloud layers are again made. Whereas the procedure developed In the 
preceding subsection (i) applied to quasi-vertical clouds such as cumulonimbus, 
the following treatment applies to multiple layers of quasi-horizontal clouds. 


O 


The Indexing scheme is shown in Fig. 11-2. Let i and j be even and 
-T&w, integers, respectively. Layer j means the layer between levels j~I and 

A 

i+1. . 

* v Define a cloud state index, 4., which can have only the values 0 and 1, 

If JU - 1 (or 0), the incident solar radiation passes through a cloudy area (or a clear 

area) in layer j. In an n-layer model, theykfafe redexUs an n-dimensional vector, 

ChucJ J} 


\ 


A 


A 
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(l. f i., I )t^ If the cloud layers are scattered, there are 2 n differed , 

possible cloud si dies. ^ 

^ °* - ^Ci /* W 

Let the optical path length of the layer between the top of the 
atmosphere, level T, and a reference level i be denoted by D . . Since 
there are t/2 layers between levels I and there are 2 ^ different possible 
cloud states for the solar radiation transferred from the top of the atmosphere 
to level u Each cloud state corresponds to a different value of D t j. Let 
D T |(V *3' •••# ^ ) denote D tI for the cloud state U x , & 3 , 

The optical path length in the layer between levels i~2 and i, 

AD. may have three different expressions corresponding to the different 

f~ 3, I , 

: cloud 'states above level ? as follows; 


(a) when solar radiation passes through a cloud area in the i-i 

(aver ~ 1), • 

' t-i 

* 

AD. . = 1 .66 u . « 

• l~3,l Cf~l / 

(b) when solar radiation does not pass through any cloud area 
in the layers above level i(4. ~ 0, 1 ^ j ^i“l), 

AD. . = (u* “ u.*) sec ^ ^ . 

|“2,I i-s i / 

(c) when solar radiation passes through a cloud area in any 
layer above level i-2 and passes through a clear area in 
layer i-1 (one or more of the 4. for 1 ^ j ^ i-3 has a 
value of 1 and l. ~ 0), 

i-l 

AD. . = 1.66 (u* - uf) 

' 1“ s , i i-2 i / 


(Hi. 34) 


(SIL35) 






Defining 


M. = 

i 


1 

for 1 <4 


(III. 37} 

3 

(i-i.) 
i=i 1 

for I 2 : 4 



A 


the cloud state for the three cases above can be classified as 

(a) / i . i = 1 

(b) &._ l = 0, M. = 1 

(c) l. = 0,_ M. 


0 


1~1 


The genera! expression for AD. . is then 


. AD. . = ]„66u . ' +(u. - u. .) [M. sec£ 

l» 3 / i Cfl-l »"1 i i 

--f K66(l-M i )](l-^ i ) , 

D t? for a cloud state (X lf . - « , 15 


T I 


* 


D T4 (°,0) = D T3 (0) + (u 3 - uj sec £ J 


* 


D T4 (0,1)-D T3 (0)+ 1,66 0^ > 


t> T4 0 #0) = D T2 0) + 1 ( u *~ u t) 


D T4 (l,]) = D f2 (l)+1.66u 


* 

C 3 


# X 


(in. 38) 


A 


(Hi. 39) ° A 



(III. 41) 


A 

A 

A 

» 

/\ 
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Let the fractional area occupied by a cloud state (X x / - • •/ ^ 

be denoted by W (Xj,*.., X^« Assuming that all cloud layers are hoi^onfa Hy^) 
quasi -homogeneous/ ~ ■ . 

' w(v.-7u^'w U) 

:;N : ~ U r ”‘ 


(111.42) 

(111.43) 


u-iy . 

where w(Z.) = (W.)(1-CL.) + X.CL. . 

I I I I 1 ; ... 

and CL. is the cloudiness of layer j. W has the property that 4,.} = ], 

/ where the summation Is taken over all possible cloud states. For example, \jvj 

fit <Z ■ 

V W(0,0) - (KL 1 )(I-CL 3 ) > 

W(0,1) = (KL 1 )CL 3 } 

- .W(1,o)'« CLj (l-CLg) ? 

W(1,1) = CL 1 *cl 3 ; . . .. 

and W(0 # 0) + W(0, 1) + W(1 f 0) + W(1 ,1) ? 1 . ' 

If the insolation incident at the top of layer [ meets a cloud area, a 
fraction R. is reflected back. It is therefore useful to define another weighting 
function W ! (x i ,..o, such that reflection is taken into account, namely 

W‘ (X, xJ = t/w*(X ) . ’(111.44) 


/ 


where 


n w (X.) ^ 


wit.) - w(X.)[(l“X.) + X.( 1 ~R.)] . 
I J I I I >* 

. =( 1 -X.)( 1 CL.)+X.CL.( 1 .-R ) 




(in. 45) 


A 

A 


A 


A 

A 


A 


A- 


A 


v\ 


The fast expression was obtained via the relation X? = (1“X.) 2 = (1-X.) ^ 

and X .O “X.) = 0. "W {^/ ? * ' ) z ^/ ) A/ - //6/^Azfe^ ^ 



Let the total albedo for both the scattered part of the radiation and 
the fraction of the absorbed part that is reflected by clouds be, for a cloud 



/ 
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state (4 R(4 lf By (111-15), with R. replaced by i.R., 1 

./ = 1 - [E(^ 1 R 1 ,. . » 7 ^g)/D(i 1 R 1 ,.o., p (HI. 46) 

Thus, for example, R(1,T) =R ia , R(0,1,0)=R 3 , R(I ,1 ,1) = R l3S and R(0,0,0) = 0. 

• The radiation budgets for an atmospheric layer and at the earth’s surface 

can now be formulated on the basis of the above results • ‘ 

. . ■ ■ i * 

i 

A bsorption In an atmospheric layer j 

The downward flux of the absorbed part of the solar radiation at levels 
t and 1+2 can be expressed as 

S .=S a £ W'u,,..., i. , )[l-A(D T ;(jt lf ..., i. ))] , (111.47) 

ai o . * 1-1 •“ i y 

m i+3 

^i+s^o S W’Ui,..*/ ^ +1 )[1“A( D T i+2 U 1 , •••/ Vi/i+i^ 

; '*•*' . • /rii jiq \ 

where the summation is taken over all the m.==2^ different possible cloud states. 
Since the net downward flux through level I of the absorbed part of the radiation 
is (W?.a.,)$ the absorption in layer i+ I is 

A,+1 al • -^CLi+i 

, (111.49} 


AS. = (1-K. )S S .. , (111.49} 

Absorption by the earth’s surface 

The absorbed part of the solar radiation that reaches the earth's surface 
is estimated by (111,47) as 

m j 1“A(D Ti (^ 1/ . * . 7 i - )) 

...... S’ T =S a 2: v X T )( 1 gTT— — — — 1 , (Hi. 50) 

aX o 1 I-i L 1-a RU,,.*-, ) J v 

s 1 

where the denominator is the correction factor due to multiple reflection between 
the clouds and the earth’s surface. However, In addition to S^, some fraction 
of the absorbed part of the solar radiation that is reflected at the tops of clouds 
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can reach fhe earth s surface after multiple reflections between cloud layers. 

Let the contribution to the downward insolation at the surface from 

radiation reflected at level i (and subsequently multiply-reflected between 

cloud layers) be * Let R y and be the total albedos of the cloud groups 

above and below level i and let the total albedo of all the cloud groups in . 

the entire atmosphere be R y + L . The absorbed part of the solar radiation 

reflected at level i Is S JtR. , 

at i+x 


After multiple reflections between the upper 
and lower cloud groups, the fraction Ru/(I-R,jR L ) has been scattered back to 
level i as a downward flux* v// Of this, the fraction (I-R L ) is transmitted through 
the lower cloud group to reach fhe earth's surface . ^Multiple reflections between 
fhe earth s surface and all the cloud groups can be introduced via the factor 
I/(I-R u+L a s ) 0 ... • 

If the i+1 layer is devoid of clouds, reflection at level i does not occur. 
Consideraing this fact, fora cloud state (j^,..., ) is given by 


CL 




S t(V 


al 
where 


RuO-Rj/C 

5 



( 111 . 51 ) 


and 


m 


( 111 . 52 ) 


I-i' ) 

. , s ;r w «v— vAA/— , 

Finally, fhe absorbed part of the solar radiation that reaches the earth's surface is 

( 111 . 53 ) 


. J %i 
S-+£ S t 
Gl i= 3 01 


4 

A 

A 

A 

A 

"a 
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Since different - values are adopted for the albedos of the scattered 

and absorbed parts of the solar radiation/ let Rj and R C^/***/ denote 

the former and R, and R(£. , . « • , l ) denote the latter . 

I 1 n 

. • 31 

In analogy to (Ul^Cf) and (UI.21)/ the scattered part of the solar 
radiation that reaches the earth's surface is 


S=$* s 1 

Si o 


i x-i ) . 

V-”'* j-AhTTit TJ > 

S C 1 i”l 


where 


a (l,-./ X T ) = l-(l-a )[1-R*U i T )3 ' .0 

C 1 . I - 1 O 41 11 

by 

'■ ' r Finally, the total absorption^? the earth's surface of both the 
scattered and absorbed pqr^of the solar radiation is 

S r ^Kl+ys'al+hJ ' (l 

The symbolic expressions developed in this section for multiple cloud 
layers appear to be relatively simple. However, the complexity and computer 
time requirement of the computation increase very rapidly with increasing 
vertical resolution. Therefore, it may be advantageous to introduce the 
following simplifications into the model. 

(a) No overlapping cloud layers. . . 

; . Overlapping of cloud layers is prohibited; at most, only one 

.. of the j&-* S is equal to one. Thus, for an n~!ayer model, there 
. ere only n + l different possible cloud states This 

■ ‘ simplification reduces the computation to a superposition or\ 

single cloud layers. ^ 
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(b) No partly cloudy condition® 

The cloudiness Is restricted to either 0 or 1 ; In any layer 
can have only a single value* Thus/ at any time and position/ 
there is only a single cloud state (i 17 . . *, 

. The radiation scheme for the current UCLA 3-Iayer GCM has been 
formulated with -fhe^o'ove--5 iTn^ T^ ^Totis-^ As shown In Appendix B, the 
scheme Is relatively straightforward 


era a), 


5 0 Surface albedo 


• ' The earth* s surface is divided Into three categories: sea surface , 

snow-free land and snow and ice surfaced, . . - . 

• . A . ; . . . 

i) Sea surface V \ . 


The albedo of She sea surface Is a function of the oyer lying cloud 
condition/ the solar zenith angle^and the wave state of the sea surface* 

. V .. Under a clear sky, a majority of the solar insolation reaches the 
earth *s surface as direct radiation. The albedo of a smooth sea surface for 
direct radiation increases from a few percent to about 100 percent as the 
solar zenith angle increases. Since the mean solar zenith angle Increases 
with latitude, the albedo of the sea surface Is usually somewhat larger in high 
latitudes than in low latitudes. . 

Under an overcast sky, almost all of the solar insolation reaches 
the surface as diffuse radiation. The albedo of the sea surface for diffuse 
radiation is nearly a constant value of 0.07* 


In the present 3-Iayer GCM, a constant value of 0.07 is adopted 

as the sea surface albedo. . 
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(ii) Snow-free land 

The albedo of snow-free land depends upon the soil condition (especially 
color) and the extent and color of natural vegetation* Coniferous forests with 
dark leaves have an albedo of about 0.14 while deserts covered by light colored 
sand have an albedo larger than 0*3* Seasonal variations in albedo caused 
by leaves withering or turning yellow are not uncommon. 

It appears that the spatial and seasonal changes in albedo are not 
essential to the simultation of the general circulation. Therefore, following 
Posey and Clapp (1964), a constant albedo of 0.14 is adopted for snow-free land. 


A 


(II f) Snow and ice 

The albedo of snow and Ice varies considerably with their cyrstalline 
structure. The latter is a function of temperature and, In the case of snow, age. 
Thus, the albedo of snow depends upon whether the snow is wet or dry, fresh or 
oidjand whether or not it is in the melting stage, the aibedo of ice depends upon 
whether^bfcW is sea packer permanent land Iceland whether or not It is in the 
melting stage* 


A 


A 


A 


In "Climates of the Polar Regions u (edited by Orvia ^1970), Putniss, 

/ * 

Vowinckel and Orvig, and Schwerdtfeger discuss^ the albedo of the snow and ice surfaces ^ 
of Greenland, the' North Polar Basin^and Antarctica, respectively. Their key results are: ^ 

(a) the albedo of snow varies from 55$ to 100$ ^ 

(b) the albedo of sea Ice varies from 20$ to 6Q$ ^ 

(c) the aibedo of permanent snow and ice fields over the interior 
of land has a lower limit of about 80$. 

In addition, the albedo of snow-covered land may vary considerably with the 
type of underlying surface. For example, the albedo of a snow-covered forest 
area depends upon the type of forest, the density of trees and the depth of the 
snow cover. 


47 



111-23 


• . . . y ' 

In midsummer, the ice and snow surfaces of the polar region receive 
more than 700 I y day of solar radiat/on. If the surface albedo changes by 
0.1, the solar radiation absorbed by a surface changes by more than 70 ly day ; 
this amount is of the same order as the net outgoing long wave radiation from 
the ‘earth's surface . Therefore, in midsummer, the radiation budget at the polar 
regions is very sensitive to the surface albedo. Since the albedos of snow and 
Ice surfaces may 'vary from 0.2 to 1.0, special attention must be paid to the 
prescription of the surface albedo for the summer polar region. 


In the current 3-Tayer model, albedos of 0.4 for sea ice and 0.7 for 
snow have been adopted, bor permanent ice and snow fields over land Interiors, 
Min(0, 85, 0.7+0.00015 h) has been employed^ where h is the land surface 
elevation in meters. 


,V "' I % VCL h 3-A ^v- &HC 




In summary j'tKef assumed values of surface albedo are: ■ 


Surface 

snow-free land 
sea • ’ 

sea ice 

temporary snow 

permanent land ice 
and snow 


Albedo, a s 

oTu ~ 

0.07 

0.4 • 

0.7 

Min(0.85, 0.7 +0.00015 h) 


/ 


M \ 

-f. i 




IV. VERTICAL DISTRIBUTION OF WATER VAPOR AND CARBON DIOX1I 


1 . Water vapor 






* The effective water vapor amount in an air column of height, z, u z / : j 
/ can be obtained from equation (II. 3) if the vertical distribution of water vapor: 

I mixing ratio q is known. In principle, since the UCLA 3-Iayer GCM predicts! 

I q for each layer from the moisture conservation equation, u* for each level 

I • L 

in the troposphere could be estimated straightforwardly from (11,3). However, ; 
j for simplicity, an assumed vertical distribution of q will be employed rather ^ 
than that predicted by the GCM 

n-be-«xprB'S5'etli5y (£/ 

— — r- ^ tTH/ ^ a p • • \y 4/ r 

yj Two methods, based upon different-formsTof this relation for the entire tropo- 
A_sphere, will be developed subsequently. 

!.V\ ' 

Y' • Although the present UCLA 3-fayer GCM does not predict the water 

vapor content in the stratosphere, this quantity is requisite, for the radiation 
* . A 

calculation. Aircraft observations by the British group (for example, Murgatroyd 

et al (1955)) have shown that the frost-point temperature decreases rapidly near 

the tropo pause and tends to asymptotically approach a constant value of about 

190°K in the lower stratosphere. Manabe and Moller (1961) summarized the 

results obtained by the British group and showed that the frost-point temperature 

in the lower stratosphere could be assumed to be 190°K, irrespective of latitude 

and season. This assumption was employed in the UCLA 2-layer moist model. 

; More recently, Williamson and Houghton (1965) estimated from radio- 
meter measurements that the stratospheric mixing ratio Is in the neighborhood of 
3x10^ g/g at least to 25 km. Mastenbrook (1968), employing balloon-borne 
frost-point hygrometers, measured the water vapor In the stratosphere over 
Trinidad, Washington, D.Cyand Thule, Greenland during 1964 and 1965. He 
found that the median vertical distribution of water vapor in the stratosphere to 


- 

".a- 


iV-2 




„ height of 28 approximates a constant^mixjng ratio within ' 

the ra nye 2. X A? 


Nil JIMAiliy IVJI'V-' 

r* 3 g/o ' 


McKinnon and Morewood (1970) also measured the water vapor distri- 
bution in the upper troposphere and lower stratosphere from 70°N to 40 S over 
North and South America by an automatic tracking solar spectrometer mounted . 
in a let aircraft. They found that the seasonal variations were within a few 
pej^cent in the lov/er stratosphere and that median values about 17.7 kryJurmg 
Northern Hemisphere winter^ gradually decrease from 1.75 x 10 at 65 N to 
1 .25 x 10” 6 at 30°$. These values arejoughiy half those obtained kyJAasfenbrook 
and Williamson and Houghton. the current UCLA 3-iayer GCM, \ 
.^coi^a itf ^ 2* 5 x I(f 6 g/g h assumed for the s^iospheve above 

r */> k J , ... 

Although relation (IV. 1) is adopted as the vertical distribution of the 
mixing ratio in the troposphere, there are many possible methods by which q 
can be estimated from (IV. 1). In the following, two methods suitable for the 
3-layer GCM wtThXe presented. . • 

(i) Method A - — 

‘ * Throughout the troposphere, the vertical distribution of q is assumed 
: io be expressed by (IV . I) with a single value of k. Thus 


^ = { 


q„ (p/p„) for P > P, T “ ?0 ° mb 

ST ST 


q ~ 2.5 x 10 for p £ p 


By (11.3), the effective water vapor amount i 






IV-3 





where p = p a+l q :Jygp^ = . 2.55^x (0.1) a x 10 ^(crn} p 

y = k+a + 1 . ’ 


r 4 Crr 


% 


& 


V 


0* 



(IV.4) 


It is now necessary to estimate a value of k at each grid point. 
From (IV.2), 

log(q/q ) =k log(p/p ) * 


‘ST 


St 


(IV. 5) 


The simplest way to obtain a single value of k for an N “layer tropospheric 
model may be ~ : « — . 

k = 1/log(q./q | >/5log(p./p ) . 

i ST , — ^ i *»T / 


r / 

= M .^(q./q ))/log( ^(Pj/p )) 

-/. r i <) f 1 J 


(IV. 6) 


/VI =. JJ/v/ 

Fora 3“Iayer model, (IV. 6) becomes 

k = log(q 1 * qa- qg/q^ViogCp/pa • ^p 
However, if a least squares method is used, S_ 

k = .^iogCq./q ) *log(p /p )/' L^log 3 (p./p ) 

r 1 J y <p k 

This expression may be superior to (IV. 6) but it also requires much more computer 
time. 



sr ? 


(IV. 7) 


(IV. 8) 


■f P : is larger than 400 mb, -the range of change of log (p./p ) is less 

than that of log(q./q ). If iog(p./p ) is replaced in (IV. 8) by some mean value, 
U st sr 

(IV* 8) reduces to'" (IV. 6) as a limiting case. 


51 



1V-4 
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(it) Method B ' \ 

Two values of k are used to express the vertical distribution of mixing 
ratio In the troposphere. That is 


q = q c (p/p c ) 1 
q = q sT ( P^T )ks 


for p.sp 

C ) 

for P st * P < P c 


(IV. 9) 
(IV.10) 


v "6 


) 


q - q =2.5x10 for p £ p = 100 mb 

... « . n n ST r r ST / 

where p Is a critical level in the middle troposphere around 400 mb. 


/ 

A 

h 


In this case, the effective wafer vapor amount can be expressed a$ 


follows: 





where 


-(W‘l for p~s 

&> r a IAp c / Vp c / J 


* P. 


> 


it it Po r / ^c \^ 3 / 

fe =U C + "(p3 ] ^Psr^Pc , (IVJ2) 

st , *•*' s 




p... 


a - a+ i / a 

Pi P c % s ^oo ^ 


(IV. 11 


' n *K +a+ i . j... 

- . .. h - p ST %/gPoo ) . - \ • 

■ r 3 = k 3 + a +1 - 

To obtain k x and q^, the least squares method is applied to the relation 
... - y. = a+k. x. 

>* 1 r > 

where x. = log(p./p ), y . = log q , a = log q , 

y y y y 


A 


A 


] fr^kr) A 




A 

A 

A 

'A 

A 

.'A 
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Then/ 


j =i J r 1 


y. )/D 

V ; 


77^y (nV SV . , v 

y q c -exp((.S y. Sx?-,E x.y. E x )/D) 

A ryr x jrtjjj = y 

/ v — ' — yijr f/\ 2 

[A ns r\y 

Using the above value of q , k 2 is obtained from 

c 

- k 2 = : !og(q c /q sT )/iog(p c /p sT ) ,, 


(IV- 3 3) 


(IV. 14) 


(IV. 15) 


The results of the three methods for estimating k, (IV. 6) and (IV. 8) 
of method Aj and (IV. 13) of method B, will now be compared to observations. 
Rodgers (1967) presented the zonal mean values of the mixing ratio in the lower 
troposphere for the northern hemisphere whlch^asmusoivcalcul a ted on the basis 
of the five year data sample of the MIT General Circulation Library. The values 
of the mixing ratio were given, at 400 , 500 , 700 , 850 and 1000 mb. The seasonal 
andJatitudirxjUrheaTi valued of k dndiho mixing ratio at 1000 mb estimated by 

y ^ ^ £ 

the three methods, and the 1000 mb observations, are in Table IV-1 . 

A 

For method B, k x and k 3 are shown by numbers without and with parentheses/ 
corresponding to values for the lower troposphere from 1000 mb to 400 mb and 
for the upper troposphere from 400 mb to 100 mb, respectively (p - 400 mb). 

; \ 

With the exception of the lowest layer in the winter high latitudes \ 
where an inversion layer predominates, the form of relation (IV. 1) is in fairly 
good agreement with the climatological data in the lower troposphere below 
400 mb. It can be considered that k a deduced by method B faithfully yields 
the actual condition for the lower troposphere. Throughout the year, k x is 
larger in midlatitudes than in the tropics and high latitudes and, with the 
exception of high latitudes, is larger in winter than in summer. 

<*5i *KtU»e kmtk 

( jeU g^c^Md */>& "> t-i 



/ 



M IV~6 

k . e ^ fa 

TABLE IV. 1 Seasonal and lafifudinalAnean valued of k in the troposphere and , 
the mixing ratio at^Me^lOOO mb/estimated by the three methods. 

• " -^B-jff|ethod B yields two values of k(k x and k 3 ). The values for ./ 

the upper troposphere between 100 and 400 mb, k 3 , are shown 
in parentheses. 


A 

A 

A 


LAT. 


k 


( 

000 mb 

(g/Kg) 


(°N) 

A 

< 

A*;- 

K&i K, 

A 

. A ! 

B 

OBS. 

•' 17 

January 

0 3.90 

' J . 

3.93 

. i 

i ! 

} ! 

2.86 (4.31) 

19.9 

21.3 

13.6 

13.8 

20 

3.65 

3.65 

3.59 (3.67) 

11.1 

n.i 

10.9 

10.4 

40 

3.26 

3.26 

3.36 (3.22) 

4.6 

4.6 

4.7 

4,4 

60 

2.85 

2.84 

2.90 (2.82) 

1.8 

1.7 

1.8 

1.5 

70 

2.65 

2.65 

2.39 (2.75) 

1.1 

1.1 

1.0 

0.8. • 

April 

0 

3.92 

3,95 

3.06 (4.26) 

20.8 

22.3 

15.2 

14.9 

20 

3.72 

3.72 

3.70 (3.73) 

13.1 

13.1 

13.0 

12.1 

40 

3.42 

3.43 

3.23 (3.50) 

6.6 

6.7 

6.2 

5.9 

60 

3.04 

3.05 

2.89 (3.10) 

2.8 

2.8 

2.6 

2.5 ^ 

70 

2.77 

2.77 

2.89 (2.72) 

1.5 

1.5 

1.5 

1.2 

July 

0 

4.04 

4.08 

2.71 (4.57) 

27.3 

30.1 

16.9 

16.8 

20 

3.98 

4.02 

2.92 (4.40) 

24.0 

26.2 

16.3 

16.0 

40 

3.81 

3.83 

3.09 (4.09) 

16.0 

16.7 

12.3 

12.0 

60 

3.59 

3.60 

3.19 (3.74) 

9.7 

10.0 

8.4 

7.5 

70 

3.46 

3.47 

3.23 (3.55) 

7.2 

7.4 

6.6 

6.0 

•October 




• ' 



0 

■ 4.02 

4.06 

2.66 (4.55) 

26.0 

28.7 

15.9 

16.2 

20 

3.90 

3.92 

3.15 (4.20) 

19.8 

20.8 

15.1 

14.6 

40 

3.54 

3.55 

3.25 (3.66) 

8.7 

8.9 

7.8 

7.5 

60 

3.17 

3.17 

3.12 (3.19) 

3.7 

3.7 

3.6 

3.2 

70 

3.01 

3.02 

2.70 (3.13) 

2.5 

2.6 

2.2 

2.0 


a 
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s 


The value of k obtained by method A shows a simpler seasonal and l 

r 6 0 

latitudinal variation. Since k Is directly proportional to the^mixlng rafio^tn 
the lower tropospher^s-4^w-pl^tke^^d4lnes^f^g'T^¥^^ k decreases 
with increasing latitude and from summer to winter. This behavior does not ( 
agree with the observations. •. v * 


In low latitudes throughout the year, and In midlatitudes in summer, 
method A overestimates k by about 1.0. As a result, the mixing ratio at 
1000 mb is more than 30 pe^j:ent larger than the observations. As an extreme 
case, the solid and dashed curves on the right side of Fig. 1V~1 demonstrate 
the difference between the vertical profiles of the mixing ratio obtained by 
method A and B, respectively, for the normal condition at the equator in July. 
The mixing ratio in the lowest 350 mb layer is overestimated, and that in the 
upper troposphere is underestimated, by method A. This may lead to an over- 
estimation of the cooling rate rate of the tropospheric column by long wave 
radiation by yielding a too small net upward flux at the earth s surface and 
a too large net outgoing flux at the top of the troposphere. 


Table IV-1 also shows the difference between k deduced by (IV. 6) 
and (IV. 8) of method A. The numbers labeled by A 1 and A refer to (IV. 6) 
and (IV, 8), respectively, fefhiAcasgw. 400 mb is the highest level for which 
there are observed values of q,^ as previously discussed, the 

difference between the A' and A values is negligibly small and (IV.6) can 

be employed in lieu of (IV. 8). , , .. _ _ 

_cr)fr£i)j >-■> 



(Jff* .. 7 As?-' ,. rr ‘ . 

Pt44f^cr^eptember^k^-7-]*/(lV.6) of. method A employed- in the- V CL 

// / s - y - y y /■> .? . /*•'' 

3-layer GCM ./ Method B Is currently undergoing testing prior to its 7-" . 

// y y - 

utilization. ! 




{/&£' 




. (Hi) Hemispheric mean normal values of mixing ratio 


/ | n order to' save computer time, it may be desirable to approximate 

' the interpolation factor ^discussed in section II .2,by constant values over the 

-v glolxd domain or throughout the year. Appropriate constant values of m can be 
sJ- obtained by substituting mean values of the mixing ratio info (11.24)- (11.27). 

For this purpose. Table 1V-2 presents the mean normal mixing ratio averaged 
over the northern hemisphere (0°-70°N) based on Rasmus's data. Since the (\ 

semi-empirical expressions for m are functions of log q, the hemispheric mean 
values shown as q represent the mean taken with respect to log q. The constants 
•that are deduced via methods A and B are also shown. 
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TABLE IV. 2. The northern ^hemispherarmecm normal values of mixing ratio ^ 
^^^ar^t^^various constants obtained by^he methods A and B <? 
The mixing ratio is shown In units of^/fecf. § » 



JAN 

APR 

JUL 

OCT 

ANN 


1000 mb 

6.0 

7.6 

13.0 - 

9.3 

8.6 . 


850 

3.9 

4.9 

8.4 

6.1 

5.6 

q 

700 

2.0 

2.6 

4.8 

3.3 

3.0 


500 

0.69 

0.86 

1.8 

1.2 

1.1 


400 

0.31 

0.38 # 

0.85 

0.53 

0.49 


k i 

3.24 

3.23 , 

2.98 

3.06 

3.14 


i 

kfj 

3.51 

3.66 

4.23 

3.94 

3.83 

B 

q 400 . 

0.32 

0.40 

0.88 

0.59 

0.51 


q -1000 

6.32 

8.05 

13.49 * 

9.70 

9.03 

A 

k 

3.44 

3.56 

3.90 

3.71 

3.65 


The annual mean vertical profile of the mixing ratio of wafer vapor in 
the troposphere as obtained from method B is 


0.51 x 10- 3 (fe) 3 ' U - 9.0& 10~ 3 {-.,?S) 3 - U P =« mb 


■{ 


■400 

2.5 x 10*A( j^-) 


~6 , d >3.83 


1000 

for 400 mb s p S: 1 00 mb 


(IV. 16) 


A 


& 


'A 


(iv) Pressure scaling factor 

To evaluate the effective wafer vapor amount, the pressure scaling 
factor a must be specified in (11.3). 
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Relies 


Elasser (1942), on the basis of theoretical and laboratory experimented^ 

^ A 

tS^udiacv adopted a value of a = 0.5 for his radiation chart. Yamamoto (1952) 

•v* I , 

recommended a ~ 1 .0 on the basis of studies presented after Elasser s paper. 

Based on the laboratory experiment by Howard, et a I (1955), Manabe and 
Moller (1 961) adopted a — 0.6. : • 

Water vapor has many absorption bands in the infra-red spectral region. 
Each band consists of many lines whose individual shapes can be important to the 
radiation transfer computation. 

Consider a single absorption line. For a weak line whose center is not 
blacked out, the total absorption does not depend on pressure broadening, hence 
a = 0. For a strong line whose center region is completely blacked out, theory 

shows that the total absorption i s proportional to the pressure, that is, a = 1 . 

• ' ■ . . 

In the actual atmosphere, absorption by water vapor occurs via an 
ensemble of weak, intermediate and strong lines^and the fractional percentage 
of each line changes with the amount of the absorber. Therefore, the value of 
a may vary from 0 to ^depending upon atmospheric conditions. 

In the lower troposphere, since p/p ~ 1 , the particular choice of 
a does not Itre^T affect the radiation budget. However, the radiation budget 
in the stratosphere may be very sensitive to the choice of a^since p/p QQ 0*1 . 
Consequently, if the pressure scaling approximation embodied in (II. 3) is to be 
utilized fora GCM which includes the stratosphere, a careful selection of a is 
■ required. . . ... 

OcL!\ c-kc, 

A value of a ~ 1 was adopted for the.2”layer moist -medoW A value of 
W i-'A UCiA ** f 

‘ a “ 0.6 is currently the^3-layer GCfA; Preliminary results indicate 

that the difference 43e4wee» the cooling rate of a tropospheric layer for these two 

0 A J 

values of Oj is about 0.1 C/day. 


A 

A 


A 


A 

A 

A 

A 
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2 • Effective amount of carbon dioxide 




!n the earth 1 s atmosphere, the spatial and seasonal variation of carbon 

<L p£<V.c*v>r* Excepir - _ p 

dioxide content Is less than 2fy<ffih Qtj gl < r T C O s i-ngra pt dly-vta— 

studies of the influence of CO^ on climateyHi*5k«*- 
‘b® assume^ that the CC^ in the atmosphere Is constant in space and time, 

it is assumedyn the followingyhat the amount of CO^ is 0.032 percent by volume 
(320 ppm), that Is 0.0489 perjspnf by weight. 1 _ 

Following Yamamoto (1952), the pressure scaling factor a for CO^ Is 
assumed to be 1.0. The effective amount of CC >2 from the earth's surface to a 
level which is usually expressed as a thickness at normal temperature (0 C) 

/. and pressure (1 atmosphere or 1000 mb), is given by dividing (11.3) by p^yp as 
! • follows: 


u co 2 ^' P s } - ’ 


SPntp t qc0 2 pfr 


dp 


V 

where 


t 


r 

(}) • 

''h' ^ jr 


= h(p* -Pp^oo ) 


q C0 2 p oo 

2g fkrp 


= 126 cm-NTP 




(IV. 17) 


(IV. 18) 


and q^Q Is the percentage content of CC^ by weight (0.0489$), P^yp <s 
^ -3 -3 

the density of CO^ at NTP(1 .977 x 10 g cm ), p^ is the surface pressure and 
.h is the total effective amount of CO^ In an air column. The effective amount 
of ££*2 in the layer between p. and p. is then ■ 


u CO = 126 |P?"Pjl cm-NTP 

^^2 1 1 J- I— 


(!V. 19) 


oo 
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/ 


If the pressure scaling factor for carbon dioxide Is a, the following general form 
is obtained instead of (IV. 19); 



• p. a+i 

,p. «’ ri 


(~\ - 

IL ~\ | 


\ p / 
r oo 

Vp / 
r oo 

6 


(IV. 20) 


A 





r 


kfp 
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V. SUMMARY 


-OAU^Ut/^ 7 ^r sAo-rf /^ 

Y<x*{U f ) +4. , he* 'b' v ty cU**r * v* // 

z/ous/y 4U ^ u j uMi&lx C-a -V} jL^ I4$f£p ">1(^->^\tY/c^L4 ! 


Simplified schemes /fe , p4^« ^ coTrlpt>taHan‘ , of f l^radit 3 Hofv-^jdgei*--ir>~f I've . i 
r f\ 

troposphere-an4H3M4re f, *e^ rfa ceph ave'beerrpT^&ni^^ 


general circulation models with //coarse vertical resolution 


<?— ■ - 


SSgmsrvm 


(i) Long wave radiation 

(a) A bulk transmission function was introduced and its value was deter- 
mined by interpolation between 1 .0 and the transmission function of the layer. 
The interpolation factors were parameterized as a function of the thickness of the 
layer, the pressure, the wafer vapor mixing ratio at fh*5"re fere nee level^and the 
lapse rates of the temperature and mixing ratio of the layer. The interpolation 
factor is linearly proportional to the thickness of the layer and varies slightly 
with the pressure and mixing ratio at the reference level. The effect of the lapse 
rates of temperature and mixing ratio is relatively small and is negligible in many 


cases. 


(b) If was assumed that any cloud occurring in a layer with a temper- 
ature less than -40°C allows transmission of a fraction of the incident long wave 
radiation. This effect was roughly approximated by reducing the cloud amount 
by the transmissivity of the cloud. 

(c) Empirical expressions for the transmission functions of v/ater vapor 

and carbon dioxide were presented and the errors of the expressions were tab- 
ulated. • . . 


A 


• (d) For simplicity, a pressure scaling approximation was adopted for the 

radiation computation- Scaling factors of either 0.6 or 1 .0 for v/ater vapor, and 
1.0 for carbon dioxide, were employed. *Rfe44mH^f^calcu!ations showed that 
the radiation budget In the troposphere was not sensitive to the choice of the value 


4 


for the scaling factor. 



/ 
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A 


A 


r\ 


y 


(ii) Solar radiation • • - 

(a) The solar spectrum was divided into scattered part 


AVI 


A 


and 4ke-ab$orbed part* The scattered part, with wavelength less than 0,9p 5 
undergoes conservative Rayleigh scatter! ng^Wh-ite absorption by water vapor 4=s*- 

The absorbed part, with wavelengths larger than 0.9 u, undergoes 

^>1 P. . . ft " 

absorption with negligible Rayleigh scattering. This € 

A 

simplifies the solar radiation calculation for a clear shy, 


, . W \G Yy r-v\ 
pt significantly 


A 


A 

A 

x 




■ (b) A cloud layer was regarded as a fictitious layer of water vapor 

that absorbs the same amount of solar insolation as the cloud layer. This con™ 

rtw 

cept considerably simplifies, the solar radiation calculation for #re~i atmosphere 

A 

with clouds- The values of the equivalent cloud water vapor amount were 
estimated for four typical cloud iypes o 

(c) General schemes for estimation of the absorption of solar 
radiation in an atmospheric (ayeijand at the earth's surface, were presented for 
a clear sky, §}> a quasi “vertical single cloud layetjand 
• *0f quasi -horizontal multiple cloud layers- . ; ' 

(■/III) Ue-Yf7^rtJ> ^ I 

“ ^ ants 

reqfftre^&yThe radiation calculation, il’^ts-weesseny"!^^ the vertical 

A /ynujf /“ 

distributions of the mixing ratios of water vapor and carbon dloxide^throughout 
the entire atmospheric column. It was assumed that the mixing ratio of water 
vapor tn the stratosphere above 100 mb has the constant value of 2-5 x 10*"^ g/g 
irrespective of season or latitude. Furthermore, it was assumed that the 

amount of carbon dioxide by volume is 320 ppm throughout the cl^bcrl- 
atmosphere. It was shown that the norma! vertical profile of the w ater va por 

mixing ratio in .the . troposphere is/fairfy, < expressed by p^cmd 

pv&cj • ^ -*r i# Cr< 

that the' vertical profile of the water vapor mixing ratio 'prsd^ted-J^ the GCM 

. ^ Aki ... A 

be adjusted to pfor the radiation calculation - 
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Appendix A: The finite difference scheme for long wave radiation for 1 

the 3-layer mode! 

The general form of the finite difference analog for the Flux of long wave 

\jLih 

radiation is given by (11.45). In the current^ 3~layer GCM, the cloud thickness 
is defined in a somewhat different manner from the general case. As shown in 
Fig. A-l , the lowest cloud is situated between levels 4 and 5.<: 

T_ =T q , we obtain / / L— 
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> -f / 

-f o-'c. , |fi ' 
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= (ctT 4 t +t AB +-22-^AB Vl-CL')(l-CLi)(l-CLi) 
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where 


A~3 

AB.. - aTt-aTt k 
M I 1 J 

■ Cl] = a.' CL. 3 

I ! I r 


and CL. is the cloudiness and a! Is the greyness of a cloud. If the temperature 
in a non-convective cloud layer is less than -40°C, a . 1 = 0,5 is assumed. The 
second term from the fasten the expression for is an additional term neces- 
sitated by the assumption that the bottom of the lowest cloud is level 5; ir^s 
denotes the interpolation factor m for the layer between levels 5 and 6 . 
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Appendix B: The finite difference scheme for solar radiation for the '■ , 

3~layer model 

In the following, the general expressions developed in section HI. 4 (ii) 
for the solar radiation calculation will be specialized to the present UCLA 
3-layer GCM* 

The cloudiness CL. is restricted to either 0 or 1 - Thus, if l, = 1 , i 

I I 

CL. - 1 and, if Z. - 0, CL. = 0; therefore, ; 

— CL. 0 (B«l) ° 

In the following, CL. will be substituted for X. everywhere. 

Employing the relations ; 

.::. V ,-V • A ■'■■■' ; ' ' ", ' '••• ‘ • 

7 (1-CL .)3) - (1-CL.) , 

( I ? A 


' - Cli? - CL. > A 

I I ' ' ' 

CL.O-CL.) = 0 j ; A 
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(i) Absorption in the atmosphere 

Combining (111*38), (111.39) and (B.3), the optical path length of the 
layer between the top of the atmosphere and level i, D t j, is: 

• (B*6) 


D to = ( u r " u 0 ) sec £ ; 


D i2 = D r 0 + 0"CL 1 )(uJ- u^)sec^ + 1.66 CLjU^ p (B.7) 

d T 4 = D IS +(i<L 9 )(u 8 *-u*)[(l-CL l )secC+l.<56CL 1 3+1.66Cl«u* #; (B.8) 

- D ts = D I4 +(I-CU.)u*[(1-CL 1 )(I-CL 3 )sece+1.66[l-0-CL 1 )(l-CL 3 )]} 

+ 1 .66 CU(u* +u ) . (B.9) 

... CS » J 

where the last term on the right hand side of (Bo 9) Ts an additional term due to 
the assumption that the base of cloud 5 is level 5 as shown in Fig. A -1. ' / 

From (III. 47) and (B.4), the downward flux at level I of the absorbed 
part of the solar radiation Is 


A 

/i 



(BUG) 

'A 

Thus, for i - 0, 2, 4 and 6, 

% 



(B.1I) 

A 

S a3 “ S o (1<R ’- )rJ ‘ A(D -- )1 9 

(8.12) ' 

A 

S a4 = s ” ( i-C | g ( i-CR 3 ) ;1- A(D T4 ) ]^ 

(B.13) 

A 

'•'VvO S a6 =S^(I-CR 1 )(1-CR 3 )(1 -CR s )[ 1-A(D ts )] # 

(B.14) 

> 

From (111.49), the absorption in the atmospheric layers \ 

= 1 , 3, 5 is 


: AS 1 =(l-CR l )S ao -S a3 , 

(B.15) 

A 

: AS a =(KR 3 )S a3 -S a4 ; 

(B.16) 

A 

A S c .= d-CR 6 ) S a 4 - S cc o 

(B.17) 

A 


70 



(il) Absorption by the earth's surface 

\ 

Let • CR. = R(CL.) = CL. R. j 

CR..=R(CL.,CL.) 

(B . 1 8) 

(B.19) 

CR... “R(CL., CL., CL.) , 

(B.20) 

Applying (111.46) to (III.1 3) ^jlfct5j) yields 


(1-CR.)(1-CR.) 

(B.21) 

{'D — I - — I- . 

^ K i« 1 “ CR.CR. • > 

■ 1 i 1 


and 


(1 — CR 4 )(1 -CR 3 )(I -cr 5 ) 

CR 13 S =1 ' l-(CR 1 CR 3 +CR a Ci^_+CI^CR^2CR l 'CR 3 CR^' 2 


= 1 "2-M-(CR 1 +CR a +CR 6 )+CR 1 CR 3 CR s > . 

where M - _(1-CR 1 )(1“CR S )0“CR S )* 

(a) The absorbed part 

• The absorbed part of the solar radiation that reaches the earth’s surface 


can be obtained from (II (.50) -(III .53). From (1 1 1. 50) and (B.20), 

.... 


s 1 = 


06 


as (l-aCR 13S ) 


(B.23) 


From (111.51), (III. 52) and (B.18)-(B.20), the contribution to the downward flux 
of solar radiation at the earth’s surface from multiple reflections due to the fop 
of cloud »+l, S* , Is 


</") \ 
It 

C/J 

a 

cr 3 

CRj (1-CR 3 s) 

(B.24) 

(l-CRjCRasKl-aCR^) ) 

S 4 - S 

as Q4 

CR S 

CRj 3 (1 rCRs ) 

(B.2 5) 

(1-CR 13 CR 5 )(1-^ s CR 135 ) / 
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The total absorbed part of the solar radiation reaching the earth's surface is then 
by (HI. 53), 


S’ +$s +S 4 * 

as as Q6 


(B.26) , */ 




(b) The scattered part 

The scattered part of the solar radiation that reaches the earth's surface 
can be obtained from (I II. 54) and (HI. 55), 

. s 1-a 


S = S — - 
$6 o. 1-a a 


c s 


: ... ... a. = Hl-a o }(l-CRj JB ) ; 


(B.27) 

(B.28) 


where CR* sS is the total albedo of the clouds in the entire atmospheric column; 
CR; 35 is found by replacing CR. by CRj in (B.22) and using CRj = CL.Rj where 
R.‘ is the albedo for the scattered part of cloud j. 


a 


A 


A 


. Finally, the total downward flux at the earth's surface is S' +S 2 + S* +S 

as as as ss 

Therefore, by (111.56), the absorption at the earth's surface is 


S 6 = (1 -a )(S 1 +S 2 + S* + $ ) , 

s as as as se ' 
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i-« s cr 135 
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CR 3 CRj (1-CR 3 s ) 


ob +S a 2 1-CRjCR 
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a* l*CR l3 CR s 
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(B.29) 


c s 


In the actual programing, clouds are indexed somewhat differently from 
the present description. Cloud types 1,2 and 3 refer to cloud layers 1,3 and 5 
of the present indexing scheme, respectively. The values of the equivalent cloud 
water vapor amount, u . , for cloud types 1, 2 and 3 are adopted for high, middle 

I 

and low clouds, respectively, as shown in Table II 1-2, The surface albedos 
currently in use are presented In section III. 5. 
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Appendix C: Required accuracy of m 

For simplicity, consider the layer between levels i-2 and I shown in 
Fig. From (II. 16) and (11. 17), the contributions of this layer to the net 

flux of long wave radiation at levels i and i-2 are respectively 


and 


C. . - (irB. -tt B. )t. . 

1-1,1 l 1“ 2 ‘~3/> ) 


. ' - (rrB. - irB. ) t . 

!“1, »“3 l 1-3' I ,1-3 <» 


From (II. 11), the fractional contribution of the above fluxes to the 
heating rate of the layer between i-2 and i, h, Is 


h = 


g(C. C. ) 
i-a,i i- 1,1-3 


c Ap 
P 


Assuming a constant temperature lapse rate, F = 31/00, h becomes 

H=4gc 1 crT? F(t. .~t. . ) 

p l 1 — 2,1 1,1-3 ) 

where the following approximation has been used 

Tf-Tt = 4T?0>T. ) = 4Ap TT? , 

- ^1 1-2 * I I 1-3 r I 

If t has a maximum error of (At the srror of h. Ah, Is 




Ah<8gc- 1 aT-r(A T ) MAX . 


The maximum possible value of Ah is then 


(Ah) 


MAX 


- 8 g c ^ a T? F (At ) 


MAX 


Considering somewhat large values of T. and I^such as T. = 300°K and 
r= 8°K/100 mb, v /e find f 


^ At) MAX • ( Ah W>/ 7,2 > n 
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where the unit of Ah is °K day ^ . 
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Now, we consider that a 10 percent error Is permissible in the 

of the heating rate# The maximum cooling rate of an atmospheric layer due to long 

wave radiation is about 2°K day ^ in a clear sky. Then the maximum permissible • 

if I o “1 

error in the heating rate,«$G:-y-be- about 0.2 K day • From (C.l), t is thus permitted 
to have a maximum error of about 0.03. 

From (11.18) or (il.19), the error of the interpolation factor m Is related 
to At as follows: — At (1+m) 3 i '[*■ ' ' 

‘ ' AT(l + m)+(l -t) ' ^ (C.2) 

The maximum permissible error in rn^for At = 0.03^ is shown in Table C~1 
.for several combinations of m and t. The combinations shown in the Table above 
the heavy solid line occur in the real atmosphere only in extreme cases when Ap 
is larger than 100 mb. From this viewpoint. It Is evident that the maximum 
permissible error in m is at least 20 percent and increases with distance below 
the heavy solid line. , _ . . 

TABLE C.l. Required accuracy In m for the accuracy of 0.03 in t. 


m 

0.1 

0.2 

0.3 . 

0.4 

0.5 

0.6 

0,7 

0.8 

0.9 

1 

0.1 

0.1 

0.2 

0.2 

0.2 

0.3 

0.3 

0.5 

0.8 

2 

0.3 

0.3 

0.3 

0.4 • 

j 0.5 

0.6 

0.7 

0.9 

1.4 

3 

0.5 

0.5 

0.6 f 

0.7 

0,8 

0.9 

1.1 

1.5 

2.2 

4 

0.7 

0.8 

0.9 

1.0 

1.2 

1.4 

1.7 

2.1 

3.0 

6 

1.3 

1.5 

1.6 - 

1.8 

2.1 

2.4 

2.9 

3.6 

4.7 

.8 

2.1 

2.3 

2.5 

2.8 

3.2 

3.6 

4.3 

5.2 

6.6 

10 

3.0 

3.2 

3.5 

3.9 

4.4 

. 5.0 

5.8 

6.9 

8.4 

15 

5.6 

6.0 

6.5 

7.1 

7.8 

8.7 

9.9 

11.3 

13.4 
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Appendix D: The seasonal variation of the sun's declination and distance from 

the earth. 

In a numerical simulation of the seasonal variation of the genera! circulation 
of the atmosphere, the seasonal variation of the solar declination, 6 , and the earth- 
sun distance, r^, must be known. 

The apparent orbital motion of the sun about the earth is shown in Fig. D-1 

by the ellipse ABA'B'A with the earth located at focus F. The circle ACA'C'A 

represents the orbit of a fictitious sun which moves with constant speed and a one 

year period. Let the sun and fictitious sun be at perigee. A, at time t and at 

o 

positions S and S' at time t, respectively. * 

The mean and true anomalies of the sun are shown by the angles M and cu. 
The mean anomaly is given by 

M W=y(t-g 

where T is the one year period. The date of perigee varies annually from Jan. 2 
to Jan. 5. The mean date of perigee for the years 1950-1972 was Jan. 3.36 (UT). 
Thus, letting t represent the time in days from the beginning of a year, t is 
2.36 days. Assuming T equal to 365 days, M is given by 

M(t) = 0.0172142 (t -2.36), radians 

An asymptotic solution of Kepler's second law in terms of the eccentricity 
e of the elliptic orbit yields 
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./ 

r ^(t)/Fg = A A ; cosM - A a cos 2M - A 3 cos 3M - . . . 


and 

a>(t) — M + BjsinM +B 3 sin 

2M +B 3 sin3M + • . 

where 

A = 1 +e 2 = 1.00027956 
o 

A 3 3 5 5 

A 1 = e - g- e " 32 e “ • * ♦ 

« 0.01671825 


A = i p 4 - 

2 e 3 e * * 

« 0.00013975 


, 3 3 135 5 

A 3 8 e 64 e ... 

^ 0.00000175 


B i = 2 e -l e3 + 4 eB " • ‘ 

^ 0.0334388 


D 5 , 11 4 , 

. B s - 4 e 24 e + 

^ 0.0003494 


g ~ 13 g 3 - -^45 e 5 + 
. 3 12 940 

& 0.00000506 


and r r is one astronomical unit. The numerical values given above are for the 
c 

earth's orbit (e = 0.01672). 

Letting £and c denote the ecliptic longitude of the sun and the inclination 
of the earth's orbit (23°27') / respectively, the solar declination 6 is given by 

6 (t) = sin" 1 (sin e sin l) 
v/here l (t) = u>(t) + 

and the ecliptic longitude at perigee, is -1 .3550737 radians (-77.64°). 6 
is positive for 0 £ X £ it and negative for ir £ l £ 2tt. 


76 


/ 



